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Abstract 



We study directed last-passage percolation on the planar square lattice whose weights 
have general distributions, or equivalently, queues in series with general service dis- 
tributions. Each row of the last-passage model has its own randomly chosen weight 
distribution. We first show the existence of the limiting time constant and list its prop- 
erties. Next we study the problem for models with Bernoulli and exponential weights, for 
which we already have more precise results. We then present some universality results 
about the limiting time constant close to the boundary of the quadrant. Close to the 
2/-axis, where the number of random distributions averaged over stays large, the limiting 
time constant takes the same universal form as in the homogeneous model. But close 
to the x-axis we see the effect of the tail of the distribution of the random environment. 
In particular we will give some estimates of the upper bound in this case. 
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Introduction 
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This paper studies the limit shapes of some last-passage percolation models in random 
environments. Specifically, we will first derive the hydrodynamic limit of the last-passage 
time for the corner growth model with exponential weights and for two Bernoulli models 
with different rules for admissible paths. Next, we will present some universality results 
for the limit shape for a broader range of underlying distributions. 

We begin by introducing the corner growth model through its queueing interpreta- 
tion. Consider service stations in series, labeled 0,1,2, ... ,£, each with unbounded wait- 
ing room and first-in first-out (FIFO) service discipline. Initially customers 0,1,2, ... ,k 
are queued up at server 0. At time t — customer begins service with server 0. Each 
customer moves through the system of servers in order, joining the queue at server j + 1 
as soon as service with server j is complete. After customer % departs server j, server j 
starts serving customer i + 1 immediately if i + 1 has been waiting in the queue, or then 
waits for customer i + 1 to arrive from station j — 1. Customers stay ordered throughout 
the process. Let X(i,j) be the service time that customer i needs at station j, and 
T(k,£) the time when customer k completes service with server t. 

Asymptotics for T{k,l) as k and I get large have been investigated a great deal in 
the past two decades. A seminal paper by Glynn-Whitt [6] studied the case of i.i.d. 
{X{i, j)}. They took advantage of the connection with directed last-passage percolation 



given by the identity 



T(k,£)=max V X(i,j) 

7T ' " 



[1.0.1) 



In this model, X(i,j) is a random weight assigned to the point The maximum 

is taken over non-decreasing nearest-neighbor lattice paths 7r C from (0,0) to (k,£) 
that are of the form 71 = {(0,0) = (x , y ), (x-i, yx ),■■■, (xk+t, Vk+i) = (k,£)} where 
(xi,yi) — = (1,0) or (0,1). Below is a picture of an admissible path from 

(0,0) to (4,3): 
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Figure 1: An admissible path to (4,3). 



It is easy to see that both the queueing setting and (1.0.1) satisfy the following 
recursive relationship for positive k and £: 



T(fc, t) = max{T(fc - 1, t),T(k, £-!)} + X{k, 



[1.0.2) 



Therefore if the process {X(i,j)} in the last-passage model has the same distribution as 
{X(i,j)} in the queueing model, then {T(k,£)} defined in (1.0.1) and in the queueing 
model have the same distribution, too. (1.0.1) together with earlier references to this 



3 

observation can be found in [B] (see Prop. 2.1). This particular last-passage model is 
also known as the corner growth model. 

Next we add a random environment to both the queueing and last-passage percolation 
models. The environment is a sequence {-F}} jeZ+ of probability distributions, generated 
by a probability measure-valued ergodic or i.i.d. process with distribution P. Given the 
sequence {Fj}j £ z + , we assume that the variables {X(i,j)} are independent and X(i,j) 
has distribution Fj. In the queueing picture this means that for each j G Z + the service 
times {X(i,j) : i G Z + } at service station j have common distribution Fj, and at the 
outset the distributions {Fj}j & z + themselves are chosen randomly according to some 
given law P. Obviously the labels "customer" and "server" are interchangeable because 
we can switch around the roles of the indices i and j. In the last-passage percolation 
model, the random environment means that weights assigned to points on the j— th row 
follow Fj. 



Although (1.0.2) is simple and clear, it does not suffice to provide much information 
about T(k,£) when k and £ are large. In fact, it is not very realistic to ask what is the 
distribution of T{k, £). Instead, we let k and £ go to infinity and scale T{k, £) in a proper 
way. The asymptotic regime we consider for T(k,£) is the hydrodynamic one where k and 
£ are both of order n and n is taken to oo. Under some moment assumptions standard 
subadditive considerations and approximations imply the existence of the deterministic 
limit for all positive real numbers x and y: 

ty(x,y) = lim n~ 1 T(\nx\, \ ny\). 



We will also verify some properties of ty(x,y) in Section 



2.2 



homogeneity, concavity 
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and continuity. 

Only in the case where the distributions Fj are exponential or geometric has it been 
possible to describe explicitly the limit ^. This is the case of ■ /M/l queues in series, 
which in terms of interacting particle systems is the same as studying either the totally 
asymmetric simple exclusion process or the zero-range process with constant jump rate. 
For i.i.d. exponential {X(i,j)} with rate 1, the limit ^(x, y) = (yfx + was first 

derived by Rost [T7] in a seminal paper on hydrodynamic limits of asymmetric exclusion 
processes. 

The random environment model with exponential F/s was studied in pfl [121 120] - 
The exact ty{x,y) can be described implicitly. It depends on the specific distribution 



of the exponential rates. In Section 3.2 we will see some explicit estimates of ^(a, 1) 
and \l/(l,a) when a is small. These two quantities have different behaviors and will be 
discussed in further details. 

Let us now set aside the queueing motivation and consider the last-passage model on 
the first quadrant of the planar integer lattice, defined by the nondecreasing lattice 
paths and the random weights {X(i,j)}. For the queueing application it is natural to 
assume the weights nonnegative, but in the general last-passage situation there is no 
reason to restrict to nonnegative weights. 

The ideal limit shape result would have some degree of universality, that is, apply to a 
broad class of distributions. Such results have been obtained only close to the boundary: 
in [13] Martin showed that in the i.i.d. case, under suitable moment hypotheses and as 
a \ 0, 

*(l,a) = /i + 2or v / a + o(v / «), (1.0.3) 
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where /x and a 2 are the common mean and variance of the weights X(i,j). The o(y/a) 
term in the statement means that lim^o or 1 ! 2 [^(a, 1) — /x — 2ay/a~\ = 0. In the i.i.d. 
case \l/ is symmetric so the same holds for \l/(a, 1). 

Our goal is to find the form Martin's result takes in the random environment setting. 
\& is no longer necessarily symmetric since the distribution of the array {X(i,j)} is not 
invariant under transposition. So we must ask the question separately for ty(l,a) and 
*(a,l). 

It turns out that for *&(a, 1), where the number of rows stays large relative to the 
number of columns, the fluctuations of the environment average out to the degree that 



our result in Theorem 4.1.1 is essentially identical to Martin's result in the homogeneous 
environment. We still have ^(a, 1) = /x + 2a^fa + o(y/a) as a \ 0, where now /x is 
the average of the "quenched" mean and a 2 is the average of the "quenched" variance. 
That is, if we let /xo = / xdF (x) and a 2 . = f(x — /x ) 2 dF Q (x) denote the mean and 
variance of the random distribution Fq, and E expectation under P, then /x = E(/xo) and 
a 2 =E(a 2 ). 

There is some evidence that we can do better than o(y/a~) for the error term. If {Fj} 
is a sequence of exponential distributions, the result 

*(a, 1) = /x + 2a^/a + 0(a) 

can be shown. For general distributions with uniform boundedness, one can achieve 

3 

o(as" e ) for any e > 0. Although not yet proved, we conjecture that 0(a) should be the 
answer even for general {Fj}. The first step should be to prove this for the homogeneous 
case. 
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The case a) does not possess a clean result such as the one above. Even though 
we are studying the deterministic limit obtained after n has been taken to infinity, we 
see an effect from the tail of the distribution of the quenched mean /xq. We illustrate 
this with the case of exponential {Fj}. Now the number not of distributions Fj is small 
compared to the number n of weights X(i,j) in each row, hence the fluctuations among 
the Fj'a become prominent. The effect comes in two forms: first, the leading term is no 
longer the averaged mean /j but the maximal mean /i*. Second, if large values among 
the row means are rare, the order of the a-dependent correction is smaller than the \/ot 
seen above and this order of magnitude depends on the tail of the distribution of /Jo- As 
an exponent characterizing this tail changes, we can see a phase transition of sorts in 
the power of a, with a logarithmic correction at the transition point. 

Intuitively, the above two phenomena suggest that when very few rows compared to 
columns are available, the optimal path makes most of its horizontal movement along 
the rows with large means very close to //*. Therefore when large means are rare, there 
are not many candidates for the optimal path, so ^(1, a) — fi* tends to be smaller. The 
other extreme is that all means are //*, i.e. they are equal. In this case we can recall 
what happens in the homogeneous case and guess the first a-dependent term may be 
xfa. We will verify this idea in the exponential model, and derive an upper bound on 
ty(l,oe) that gives the correction of order yfa as well for general distributions under 
sufficient conditions. 

The key idea in proving universality results in this paper is to compare limiting time 
constant in models with general distributions to that in models with normal distribu- 
tions. For this purpose we need to quantify the difference between *$!f(%, y) and *$?g(%, y) 
for two processes {ij-}jez+ an d {Gj}jez + - An example of this is Lemma El In the proof 



we use as auxiliary results bounds on the limits of last-passage models with Bernoulli 
weights . 

It is worth noting that with Bernoulli weights the limiting time constant ty(x,y) 
has not been derived for the standard corner growth model. ty(x,y) can be solved in a 
model with Bernoulli weights when the path geometry is altered suitably. The model 
we take up is the one where the paths are weakly increasing in one coordinate but 
strictly in the other. There are two cases, depending on which coordinate is required to 
increase strictly. If we require the x-coordinate to increase strictly then an admissible 
path {(x ,y ), (zi,yi),.. . , (x m ,y m )} satisfies 



x i+ i - Xi = 1 and y < y ± < ■ ■ ■ < y r , 



(1.0.4) 



We give a figure below showing a possible path: 




2 3 4 5 
Figure 2: An admissible path to (5,4) 
with x-coordinate strictly increasing. 



The other case interchanges x and y. These cases have to be addressed separately 
because the random environment attached to rows makes the model asymmetric, i.e. 
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the value of ^(x, y) changes when we interchange the two coordinates. The sum of these 
two last-passage values gives a bound for the case where neither coordinate is required 
to increase strictly in each step. 

We derive the exact limit constants for Bernoulli models with both types of 
"strict /weak" paths. For one of them this has been done before by Gravner, Tracy and 
Widom [9]. Their proof utilizes the fact that the distribution of T(k,£) is a symmetric 
function of the environment in the sense that it is not affected if we interchange the 
distributions in any two rows (at least for the particular Bernoulli case they study). The 
proof here is completely different. It is based on the idea in [19J where the limit for the 
homogeneous case was derived: the last-passage model is coupled with a particle system 
whose invariant distributions can be written down explicitly, and then through some 
convex analysis the speed of a tagged particle yields the explicit limit of the last-passage 
model. This same approach can be adapted to the random environment case so that 
results in [19] can be generalized. 

Further remarks on the literature. In p], a different asymptotic regime given by 
-^+ sr T(\n a x\ , n) with < a < 1 was studied in the homogeneous model. They derived 

n 2 

an asymptotic result for the above quantity when the underlying distribution has an 
exponentially decaying tail. It would be interesting to see whether the result can be 
generalized to the random environment case. 

Many papers also addressed questions of fluctuations. For the last-passage model 
with i.i.d. exponential or geometric weights, the distributional limit with fluctuations 
of order n 1//3 and limit given by the Tracy- Widom GUE distribution was proved by 
Johansson [10J. As for the shape, universality has been achieved only close to the 
boundary, by Baik-Suidan [2] and Bodineau-Martin [3]. 
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Fluctuations of the Bernoulli model with strict/weak paths and homogeneous weights 
were derived first in [TT] and later also in [7j. For the model in a random environment 
fluctuation limits appear in [91 [8] . 

On the lattice Zl we can imagine three types of nondecreasing paths: (i) weak-weak: 



both coordinates required to increase weakly, the type used in (1.0.1); (ii) strict-weak: 



one coordinate increases strictly, as above in (1.0.4); and (iii) strict-strict: both coor- 



dinates increase strictly so an admissible path {(x ,yo), (xi,y%), . . . , (x m ,y m )} satisfies 
xo < ■ ■ ■ < x m and yo < • • • < y m . As mentioned, with Bernoulli weights the strict-weak 
case is solvable but the weak-weak case appears harder. The third case, strict-strict, is 
also solvable with Bernoulli weights. The shape was derived in [TS] and recent work on 
this model appears in [5]. 

Organization of the paper. We begin by introducing the last-passage time percolation 
model in Chapter [2] and verify the existence of the limiting time constant in Section 



2.2 



Next we present some results specifically for Bernoulli models (Section 3.1) and 



exponential models (Section 3.2). Then we will show universality theorems on the shape 



close to the boundary in Chapter El in Section 4.1 we present Theorem 4.1.1 on ^(a, 1 



and in Section 4.2 we have some estimates on \l/(l, a). 



Some frequently used notation. We write 



ess sup / = inf{s G R : P(/ > s) = 0} 

v 

for the essential supremum of a function / under a measure P. Z + = {0,1,2,...}, 
N = {1,2,3,...}. 1(A) is the indicator function of event A. 
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Chapter 2 

The existence and properties of the 
limiting time constant 

2.1 The last-passage percolation model 

We give a precise definition of the last-passage model in a random environment. Let P 
be a stationary, ergodic probability measure on the space .M of sequences of Borel 
probability distributions on R. E denotes expectation under P. For some of the results 
in the following chapters P will be further assumed to be an i.i.d. product measure. A 
realization of the distribution- valued process under P is denoted by {Fj}j e z + . This is 
the environment. Given {Fj}, the weights {X(z) : z G Z^} are independent real-valued 
random variables with marginal distributions X(i,j) ~ Fj for G Z+. Let (Q, T, P) 
be the probability space on which all variables {Fj, X(i, j)} are defined, and denote 
expectation under P by E. 

A (weakly) nondecreasing path is a sequence of points z = (xo, yo), z\ = (x±, y\), . . . , 
z m = {xm,ym) in Z 2 + that satisfy x < x x < ■ ■ ■ < x m , y < yi < ■ ■ ■ < y m , and 
\x i+ i — Xi\ + \yi+i — yi\ = 1. For z±,z 2 G with z\ < z 2 (coordinatewise ordering), 
let H[zi, Z2) be the set of nondecreasing paths from z\ to z 2 - Whether the endpoints z\ 
and z 2 are included in the path makes no difference to the limit results below, but to 
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be precise let us include z 1 and exclude z 2 , so that we can run a subadditive argument 
later. 

Remark 1. We included the endpoint in the definition above Figure^l\ because we wanted 



it to be consistent with (1.0.2). 

Hereafter we will always exclude the endpoint when talking about a path between two 
points. 

The last-passage time T(zx, z 2 ) from z\ to z 2 is defined by 

T(zi,z 2 )= max > X(z). 

nen[zi,za) 

2£7T 

When zi = Q abbreviate U(z) = n[0, z) and T(z) = T(0, z). 

T(z) is a random variable that depends on the underlying distributions, and will 
be quite complicated as z moves far away from the origin. However, the following 
quantity, known as the limiting time constant, exists under proper conditions and provide 
information about the last-passage time 



ty(x,y) = lim — T([nx\, \ ny\). 

n— >oo n 



2.2 The existence and properties 

We now give a set of sufficient conditions for the aforementioned limit to exist. Put 
these three assumptions on the model: 



E\X(z)\ < oo, 



(2.2.1) 



f°° r i 1 / 2 

J |l-E(F (x))| dx < oo, 



and 



PCX) 

/ esssup(l — F (x)) dx < oo. 
Jo P 

We start with these assumptions and consider the existence of ^(x,y). 
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(2.2.2) 
(2.2.3) 



Proposition 2.1. Assume P is ergodic and satisfies (2.2.1), (2.2.2) and (2.2.3). Then 
for all (x,y) G (0, oo) 2 the last-passage time constant 



ty(x,y) = lim — T([nx\, \ ny\) 

n— >oo fl 



(2.2.4) 



exists as a limit both P -almost surely and in L 1 (P). Furthermore, ty(x,y) is a homoge- 
neous, concave and continuous function on (0, oo) 2 . 



Assumption (2.2.2) is also used for the constant distribution case, see (2.5) in 



Some further control along the lines of assumption (2.2.3) is required for our case. For 



example, suppose 1 — Fj(x) = e~^ jX for random ^ e (0,oo). Then (2.2.3) holds iff 
ess infp(£ ) > 0. If the distribution of £o is n °t bounded away from zero, n~ l T(n, n) — > oo 
because we can simply collect all the weights from the row with minimal £j among 



{£q, • • • , £ n }- However, assumption (2.2.2) can be satisfied without bounding £ away 
from zero. 

Proof. We first prove the theorem for integer pairs (x,y), and then extend to rational 
numbers and finally real numbers. 

Step 1: consider (x,y) G N 2 . Set Z m ^ n = —T((mx,my),(nx,ny)) for < m < 
n, and verify that under the distribution P, Z m ^ n satisfies assumptions (i), (ii) and 
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(iii) in Liggett's version of the subadditive ergodic theorem [3J p. 358]. In particular, 
Zo, m + Z m ^ n > Z 0jn , {Z n k,( n +i)k,n > 1} is ergodic for each k, and the distribution of 
{Z m ,m+k, k > 1} does not depend on m. 

We need to work harder on condition (iv), i.e. we need to show EZ^ < oo and for 
each n, EZ 0jn > for some 7 > — 00. It is easy to see 



EZ+ < E|T((0,0), (x,y))\ < E ^ \X(i,j)\ = (x + l)(y + 1)E|X(0,0)| < 00. 

0<i<x,0<j<y 



Next we show EZ 0i „ > 'jn for some 7 > —00 under (2.2.2) and (2.2.3). This is triv- 



ially true for a Bernoulli model where given {Fj} the weights have marginal distributions 

P(X(i,j) = 1) = 1 - Fj(u) = 1 - P(X(i,j) = 0). (2.2.5) 

Therefore this Bernoulli model satisfies all conditions of the Subadditive ergodic theorem, 
and ^ Ber[i-F(u)](x, y), the limiting time constant, is well-defined P — a.s. and in -^ 1 (P) 



for (x,y) G N . We will see an upper bound (3.1.8) for ty(x,y) of the Bernoulli model 
in the next chapter. We use it here without proof in the following calculation: 

11 1 f°° 

-EZ 0n > E max Vx(z) + = E max V/ I(X(z)>u)du 

n n ireU(nx,ny) ^— ' n iren(nx,ny) ^-^ J n 

z£tt z£ir 

1 r°° 

> E / max yl(X(z) > u) du 

n Jq nell(nx,ny) 

/*oo j /*oo 

= — / sup — E max H-^( z ) > u ) ^ u = ~ / ^Ber[i~F(u)](x,y) du 

J n U Trell(nx,ny) J Q 
/*oo /»oo 

>-(y + 4y/xy) I y/l -EF (u) du - x / (l - essMF {u))du. 
Jo Jo p 

(2.2.6) 
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Here 1(A) is the indicator function of event A. By assumptions (2.2.2) and (2.2.3), 
EZo jn > n7 for a constant 7 > —00. These estimates justify the application of the 
subadditive ergodic theorem. So now for (x, y) G N 2 , we can define the following P — a.s. 
and L X (P) limit 

*ff(x,y) = lim —T(nx,ny). 

n-»oo fl 

Step 2: take (x,y) G (Q fl (0, oo)) 2 . Let x = ^ and y = ^ be in their reduced 
forms, i.e. X{ and yi are positive for i = 1, 2 and gcd(xi,X2) = gcd(yi,y2) = 1- Let be 
the least common multiple of X2 and y2, so (fcr, %) G M 2 . 

For every positive integer n, write n = Mk + r for integers M and r such that 
< r < — 1. Then we have 

Mkx < [nx\ < (M + l)kx, Mky < [ny\ < (M + 1)%. 

So if we denote zi(n) = (Mkx, Mky) and Z2(n) = ((M + l)kx, (M + l)ky), we have 
the following inequalities from superadditivity: 

T(z 1 (n))+T((z l (n),([nx\, [ny\)) < T([nx\, [ny\) < T(z 2 (n)) -T(([nx\, [ny\),z 2 (n)) 
Obviously, 



T(zi(n), ( \nx\ , [ny\ )) > — max y \X(z 



n£n[zi(n),Z2(n)) 

and this leads to 



1 j 



T{zi(n))- max V|I( 2 )|<T(H,W). (2.2.7) 

well[zi(n),Z2(n)) * — ' 



Z&TT 



Let £ > be any small positive number, 
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VP( max S^\X(z)\>ne) 
^— ' Ven[ 2l ( n ), Z2 (n))^' 

n=l zGti" 



= VP( max y^\X(z)\>ne) 

^ \ell[0,(kx,ky)) ^ 
n=l z£tt 

i r°° 

<- P( max S^\X(z)\ > x)dx (2.2.8) 
e J \enio,(kx,k y )) ^ 

= -e( max VlX^)^ 
<-(fcx + 1)(% + 1)E|X(0, 0)| < oo. 

By Borel-Cantelli Lemma, - max 7rg n[ 21 („), 22 („)) Ylzeir l^( z )l —'•OP — a.s. Dividing 



through by n and taking limit in (2.2.7) gives that 



—ty(kx,ky) = lim —TizAn)) < liminf — T( I nx\ , InyJ) P — a.s. 

k n->oo n n->oo n 



Similarly, we can show the other direction 



lim sup — T(\nx\, \ny\) < lim — Tizzirif) — — ^(kx, ky) P — a.s. 

ra->oc 71 n->oo n k 



This shows that 



lim — T( , L^J/J ) — -j^{kx, ky) P — a.s. 



The definition of has now been extended from integer points to rational points 
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by 

*(x,y) = ^{kx,ky), (2-2.9) 
where k is defined at the beginning of Step 2. 



From (2.2.7), we get 

( — TizAn)) — —ty(kx, ky)] max \X(z)\ < —T(\nx\,\ny\) — —ty(kx,ky). 

z£tt 

Similarly, 

-T([nx\, \ny\)-\^(kx,ky) < (-T(z 2 (n))- T ^(kx,ky)) + - max y^\X(z)\. 
n k \n k J n 7renUi(n),« 2 (n)) 



Note that for i = 1,2, because (kx, ky) E N 



lim E\-T(zi(n)) - \^(kx, ky)\ = 0. 

ra-»oo n k 



In addition, 



— e( max 

n \7ren[* 1 (n),2 2 (n)) 



,z 2 (n)) « — ' / n 



which also converges to as n goes to infinity. 
Therefore the -^ 1 (P) convergence follows 



lim E| -T( [nx\ , [ny\ ) - T ^f(kx, ky) \ = 0. 

n-5>oo n k 



One can fairly easily check the following properties: for positive rational pairs (x\, yi) 
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and (x 2 ,y 2 ) 

1. homogeneity: ty(cxi,cyi) = c\&(xi,yi) for any positive rational number c. 

2. superadditivity: V(xi +x 2 ,yi + y 2 ) > #(zi,yi) + #(2:2,2/2) • 

3. The above two together imply concavity: for rational < c < 1 and let x = 
cxi + (1 - c)x 2 , y — cyi + (1 - c)y 2 , we have 

*(x,y) > c^(x 1 ,y 1 ) + (1 - c)*(x 2 ,y 2 ). (2.2.10) 

Step 3: we extend the definition to (x,y) G (0,oo) 2 . First, we prove lim n T( - n '^ iy ^ 
exists P — a.s. for all y G (0, oo) 

If y is not rational, we can pick y±,y 2 G Q D (0, oo) such that y\ < y < y 2 . By picking 
the optimal path from the origin to (n, |_ n 2/iJ) an d then moving directly to (n, [ny\), we 
have the following inequality: 

\ny\-l 

T(n,[n yi \)+ J2 X(n,j)<T(n,[ny\). (2.2.11) 

j=[ny 1 ] 

We now take a random subsequence {rik, k — 1, 2, ...} such that 
— T(n fc , [rajfcJ/J) ->■ liminf -T(n, |ra/J). 

By the strong law of large numbers \ ^™2/J-l™2/iJ X(0, j) converges to (y — y\)EX(0, 0) 
P — a.s. and in probability. Therefore if we fix an e > 0, then for every £ = 1,2,... we 
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can find an integer N(£) such that 



\ny\-\nyi\ 

P (l^ E ^(0,j)-(j/-J/i)EX(0,0)|>e)<2- < (2.2.12) 

for all n > N(£). 

Since I Y}f=\_ny\\ X ( n >j) has the same distribution as i £L"i/J-L«wJ X(0, j), ( gg32] ) 
implies that from {r^} we can select a further subsequence {n^)} such that 

oo . L n fc(i)fJ _1 

E P (|— E *(r^)-(3/-3/i)EX(0,0)| >e) <oo, 

' =1 i=L n fe(()2/lJ 

and this shows that P — a.s. we have 

lim V X(n,j) = (2/-2/i)EX(0,0). 

j=l n k(i)yii 

Hence by dividing through by n and taking limits along this subsequence {n*^)} in 



(2.2.11) we get 



+ (y-yOEX^O) < liminf -T(n, LnyJ). (2.2.13) 

n 



Similarly, we can show 



limsup-T(n, [ny\) < j/ 2 ) - (y 2 - y)EX(0, 0). (2.2.14) 

n 



Note that both of the above inequalities hold P — a.s. 

So now it is natural that we want to let y\ and y 2 approach y from both sides. We 



19 

need the following lemma. 

Lemma 2. If f(x) is a concave function defined on Qfl (0, oo), then it can be extended 
uniquely to a continuous function on (0, oo) by 



f(x)= lim f(y). (2.2.15) 



Proof. Let x > be a fixed real number. We first prove the one-sided limit \im y£ Q ty ^ x - f(y) 
is well-defined. If this is not the case, then we can find two sequences of rational 
numbers {u n } and {v n } such that they both approach x from below, and the limits 
A\ = lim n /(tt n ) and A 2 = lim n /(f ra ) both exist with A\ > A 2 . Take e > small 
enough. We can find three rational numbers u, v! G {«„} and v G {v n } such that 
u < v < u', and f(u), f{u') > A\ — e, f{y) < A 2 + e. Take a rational number < q < 1 
such that v — qu + (1 — g)w' , then 

/(«) < A 2 + e < A, - e < qf(u) + (1 - q)f(u'). 



This contradicts (2.2.10), so the left limit A = \im y( zQ y ^ x - f(y) exists. Similarly, 



B = lirn^gQ y -+ x + f(y) also exists. We only need to show A = B. 



Assume A > B, and choose < e < ^y^. Then there exists 5 > such that for any 
rational numbers u and v with x — 5 < u < x and x < v < x + 5, 

A-e < f(u) < A + e, B - e < f(v) < B + e. 



Take rational numbers u and v such that x—5 < u < x < v < x+5 and ^-^ < 1 — -^%, 
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and pick a rational number a G (fz^, 1 — ~£zb)- It follows that 

af(v) + (1 - a)/(«) > a(£ - e) + (1 - a) (A - e) = aB + (1 - a) A - £ > B + e. 
However, by the choice of a, av + (1 — a)w is a rational number in (x,v), so 
/(aw + (1 - a)u) < B + s < af(v) + (1 - a)f(u). 



This again violates (2.2.10), and by contradiction we reject A > B. Similarly, we 
can show A < B is not possible either. Hence A = B and \im y£ Q^ x f(y) exists. 

For x G Q PI (0, oo), we need to show this limit is consistent with the original value 
f(x). We can repeat the above proof by contradiction and modify it when necessary. 
Specifically, we let A = f(x), B = lim 2/6 Q jJ/ _ >a .+ f(y) and assume A > B; in the following 
part we choose u = x and v G (x, x+5) such that f(v) G (B—e, B+s), and take a rational 
number a G (0, 1 — -sttb)- We can check /(av + (1 — < af(v) + (1 — a)f(u), so it 



contradicts (2.2.10) and A > B is rejected. Similarly we reject A < B and get A = B. 



Since the two-sided limit lim^gQ^a; /(y) exists, this gives (2.2.15) for x G Q fl (0, oo). 

One can quickly check that the extension keeps concavity: for x and y in (0, oo) and 
< c < 1, take rational sequences x n — > x, y n — > y and c n — > c, then 

f(cf(x) + (1 - c)f(y)) = hm f(c n f(x n ) + (1 - c n )f(y n )) 

> limc„/(x n ) + lim(l - c n )f(y n ) 

n n 

= cf(x) + (l-c)f(y). 
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Since a finite concave function on an open set is continuous by Theorem 10.1 of [15J, we 
get continuity. The uniqueness is trivial because any continuous function must satisfy 

□ 



(2.2.15) 



Let us return to the proof of Proposition 2.1 Since y) is concave function defined 



on Q n (0, oo), it can be extended to (0, oo) by (2.2.15). We let y\ and yi approach y in 



(2.2.13) and (2.2.14), and get 



lim —T(n,\ny\) — lim ^(1, u) P — a.s. 



Therefore we can extend the definition to y e (0, oo) 



= lim 



(2.2.16) 



We then extend the definition of *$f(x,y) to any (x,y) G (0, oo) 2 , and show that 



lim — T( \nx\ , I ny\ ) = x^(l, — ) P — a.s. 

n->oo n x 



(2.2.17) 



If we write m = \nx\ , then m < nx < m + 1, hence |_ m f J — \_ n y\ — L( m + 1)~J- ^ 
is clear that 



[nyj-l 



T(m, L«^— J) — |X(m, i)| < T([nx\, \ny\) 

*=L m |J 

L(m+l)*J-l 

<T(m, L(ra + 1)-J)+ ^ 

i=[nyj 



(2.2.18) 



Similarly to (2.2.8), we can run a Borel-Cantelli argument and claim that as n goes 
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to infinity, 



1 L-yJ 1 L(m+i)*J 

- V \X(m,i)\-tO and - V \X(m,i)\ P - a.s 
n n z — ' 



Therefore dividing through by n and taking limits in (2.2.18) gives that P — a.s 



— ) < liminf — T( \nx\ , \ny\ ) < limsup — T( [rixj , [_^2/J ) < x\E'(l, — ] 



So (2.2.17) is proved. We can then define for (x, y) G (0, oo) 2 that 



ty{x,y) = x^(l, 



x 



(2.2.19) 



Finally, we prove -^ 1 (P) convergence for (x,y) G (0, oo) 2 . From (2.2.11) and its 
counterpart in the other direction, we get 



[ny\-l 

(-T(n,Lnj/iJ)-*(l,yi)) + (*(l,!/i)-*(l,J/))+- £ X ( n >^') 

J=L"2/iJ 



<-T(n, |nyj)-tt(l,y) 



L'/ty2j — l 

<(-r(n,Ln3/ 2 J)-*(l,y 2 )) + (*(l,2/ 2 )-*(l s 3/))+- ^ X{n,j) 



We have shown that for rational numbers y\ and y^ lim^oo E| -T(n, \nyi\)—^f(l, y^ 
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0. Then 



lim E\-T(n, \ny\) - #(1,7/)| 

n— >oo n 

<(|*(l,j/ 2 ) - + 0/2 - 2/)E|X(0,0) 

V - + (t/ - yi)E|X(0, 



We let ?/i and t/ 2 approach 7/ and get 



lim E\-T(n, \ny\) - #(1,2/) I = 0. 



We can use a very similar argument starting from (2.2.18) to get 



lim El— T(\nx\, \ny\) — x^(l, — ' 



0. 



So iv 1 (P) convergence is proved. 

So now we have started from the definition of *f?(x, y) on integer points and extended 



it to (0, oo) 2 by (2.2.9), (2.2.16) and (2.2.19). We can immediately extend the homo- 
geneity, superadditivity, and concavity conditions to real points. Again by Theorem 10.1 
of [15] a finite concave function on an open set is continuous, we get continuity. 

□ 



Now we have finished the proof of Proposition 2.1 



We may also define ty(x, y) using supremum. If we denote x n = \nx\ and y n = [ny\. 
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by superadditivity, we have 

T{x mi Dm) + T{ix m , 2/ m ), {x m + x n , y m + 

We note that x m+n — x m — x n — or 1, and so is y m+n — y m — y n - If EX(0, 0) > 0, 
then we can easily check that ET(0, 1), ET(1,0) and ET(1, 1) are nonnegative. This 
gives ET(x m ,y m ) + ET(x n ,y n ) < ET(x m+n ,y m+n ), which leads to another definition 
that 



fy{x,y) = lim -ET(x n ,y n ) = sup -ET(x n ,y n ) = sup -ET([nx\, [ny\). 

n-^oo n n n n n 

This alternative definition will be helpful in some settings where we need to estimate 
the upper bound of last-passage times. It may not be true if EX(0,0) < 0. An easy 
counterexample is the case where X(z) = — 1 for all z £ 71?+. We easily see that 
^~ET([nx\, [ny\) = -(—[nx\ — [ny\) > ^(—nx — ny) = —x — y = ty(x,y) when x and 
y are not integers. However, if x and y are both integers, ty(x,y) = sup n ^ET(nx,ra/) 
is a valid definition regardless of the sign of EX(0, 0). 
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Chapter 3 



Results for Bernoulli and 



exponential models 



3.1 Bernoulli models with strict- weak paths in a 
random environment 



As we have seen in the proof of Proposition 2.1, last-passage models with Bernoulli- 
distributed weights can play an important role when we study general models. As for 
models with Bernoulli weights, one of the major difficulties is that there is no explicit 
results so far about *f>(x,y) with the weakly increasing paths. For this reason in this 
chapter we first study Bernoulli models with two different types of admissible paths, and 
eventually give an estimate of ty{x,y) with the weakly increasing paths. 

The environment is now an i.i.d. sequence {pj}j G i + of numbers pj G [0, 1], with distri- 
bution P. Given {pj}, the weights {X(i,j)} are independent with marginal distributions 
P(X(i,j) = 1) = pj = 1 — P(X(i,j) = 0). We consider two last-passage times that 



2 



differ by the type of admissible path: for z^^gZ 



T^(zi, Z2) = max N X(z) and T+(zi, z 2 ) = max N X(z). (3.1.1) 

7reII-).[zi,Z2) / — ' 7ren t [zi,z 2 ) 

ZS7T ZGTT 
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In terms of coordinates denote the endpoints by zj, = (a,k,bk), k = 1,2. Admissible 
paths it G IL-Jzi, z-i) are of the form tx = {(ai, yo)(ai + 1, . . . , (a 2 — 1, y a2 _ ai _i)} with 
b\ < < yi < ■ ■ • < y a2 -a 1 -i < &2- Please see Figure [2] Again note that now we always 
exclude the end point as was declared in Remark [TJ 

Symmetrically paths 7r G Ilf [zi, z 2 ) are of the form 7r = {(xq, b±), (xi, b\ + 1) . . . , 
(a;6 2 _6 1 _i, 62 — 1)} with a± < x < xi < ■ ■ ■ < x^-fri-i ^ a 2- Thus paths in II_^[^i, ^ 2 ) 
increase strictly in the x-direction while those in n^!, z 2 ) increase strictly in the y- 
direction. 

As before we simplify notation with T^(0,z) = T^(z). The almost sure limits are 
denoted by 



ty-+(x,y) = lim — T_^( [nx\ , [ny\ ) and ty+(x,y) = lim —T^(\nx\,\ny\) (3.1.2) 

n— >oo n ' n— ¥00 TL 



for (x,y) G (0,oo) 2 . The proof of the existence of the above limits can be outlined as 
follows: we first verify the assumptions of the subadditive ergodic theorem for (x, y) G 
N 2 , and note that the moment assumption is trivial for Bernoulli models because of the 
uniform boundedness; then we extend the definition to all (x,y) G (0, oo) 2 using the 
same argument as we had in the previous chapter. We will omit the details and claim 
the existence of the limits. 



Remark 3. In (2.2.6) we used a result (3.1.8) from this chapter. To remove the concern 



for circularity, we note that the proof of Proposition 2. 1 works for the Bernoulli models 



even without knowing (3.1.8). The logic progression actually should be: Proposition 2.1 



holds for Bernoulli models, then we derive (3.1.8), and apply (3.1.8) to prove Proposition 



2. 1 for more general models under moment conditions. 
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The next theorem gives the explicit limits. (3.1.3) is the same as in [HI Thm. 1]. Inside 

the E[ - • •] expectations below p is the random Bernoulli probability. Let b = ess sup p 

p 

denote the maximal probability. 

We prove the formulas and inequalities first for and then for ty^. It is convenient 
to assume b < 1. Results for the case 6=1 follow by taking a limit. 



Theorem 3.1.1. The limits in (3.1.2) are as follows for x,y G (0, oo 



bx + y{l-b)E , x/y>E 



p(i-p) 

(b- P ) 2 



1=2 

(zo-p) 2 



E 



p 

1-p 



< x/y < E 



p(i-p) 

(b-P) 2 



< x/y < E 



i-p 



(3.1.3) 



with zq G (6, 1) uniquely defined by the equation 



x/y = E 



pQ- -p) 

(z - p) 2 



y t (x,y) = < 



y - yz 2 E 



(zo +p) 2 



< x/y < E 
x/y > E 



i-p 
p 



i-p 



mt/i G (0, oo) uniquely defined by the equation 



(3.1.4) 



x/y = E 



Pi} -P) 

(zo+p) 2 




(b) Approximation of the function ^(a, 1) from simulation. 
Figure 3: simulations to verify Theorem 3.1.1 
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To illustrate this result, we present two plots based on simulations: the rates pj are 
chosen to be i.i.d. with F(pj < x) = 1 — (^|^) 3 supported on [0.4,0.9]. 

The dashed curves are the precise values described by Theorem |3.1.1 The solid 
curves are the approximations ^ ^ na ^ - and Tt ^" tt ^'^ respectively when n = 20, 000. 



In Figure 3.3(a), the approximation is not very accurate when a is close to after 



zooming in. The reason is that the density of Bernoulli rates near b is low, so the first 



few rows are not likely to have means close to b. In Figure 3.3(b), the error is almost 
negligible. We can also clearly see that ^5/^(1, a) and ^(a, 1) approach to different 
limits as a \ 0. 

Our second result gives simplified bounds that are useful for the proof of the Theorem 



4.1.1 in the next chapter. Let p = E(p) be the mean of the environment. ty(x,y) is the 



limiting time constant with weakly increasing paths defined in Proposition |2.1| 
Theorem 3.1.2. The following three inequalities hold for the Bernoulli model: 

y) < bx + 2^/p(l -b)xy, (3.1.5) 

^^(x, y) <py + 2a/p(1 - p)xy (3.1.6) 

and 

ty(x, y) < py + 4a/p(1 — p)xy + bx. (3.1.7) 



30 

(3.1.7) follows from (3.1.5) and (3.1.6) because *&(x, y) < ^/^(x,y) + ^^(x,y). An- 
other loose estimate we will use later following (3.1.7) is 



^(x, y) < py + 4a/p(1 — p)xy + bx < {y + 4 v /xy) A /p + bx. 



(3.1.8) 



Proof of (3.1.3) and (3.1.5). We adapt the proof from [19] to the random environment 



situation and sketch the main points. Results from [19j are directly applicable in the 
random environment and will be quoted without further explanation. 

Consider now the environment {pj} fixed, but the weights X(i,j) random. For 
integers < s < t and a, k define an inverse to the last-passage time as 



r((a, s), k, t) = min{Z e Z + : T_*((a + 1, s + 1), (a + I + 1, t)) > k}. 

For the special case k = we define T((a, s),0,t) = 0, but T((a,s),k,t) > for k > 0. 
Knowing the limits of the variables T is the same as knowing By the homogeneity 
of it is enough to find h(x) = ^^(x, 1). By the homogeneity and superadditivity of 
h is concave and nondecreasing. Let g be the inverse function of h on M + . Then g 
is convex and nondecreasing, and (7.4) in [19J still holds here: 

tg(x/t) = lim -r((0,0), [nx\, [nt\). 

To find these functions we construct an exclusion-type process z(t) = {^(t) : k e Z} 
of labeled, ordered particles Zk{t) < Zk+i{t) that jump leftward on the lattice Z, in 
discrete time t E Z + . Given an initial configuration {^(0)} that satisfies Zi_i(0) < 
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Zi(Q) — 1 and liminf i _ s ._ 00 |i| _1 Zj(0) > —1/6, the evolution is denned by 

z k (t) = inf {zj(0) + r((zi(0), 0), k - i,t)}, keZ,teN. (3.1.9) 

i:i<k 

It can be checked that z(t) is a well-defined Markov process, in particular that z k (t) > 
— oo almost surely. These claims are identical to Lemma 5.2 and Lemma 5.3 in [19] . 

Define the process {i]i(t)} of interparticle distances by rji(t) = z i+ i(t) — Zi(t) for i G Z 
and t G Z + . By Prop. 1 in [12] process {r)i(t)} has a family of i.i.d. geometric invariant 
distributions indexed by the mean u G [1, b~ l ) and defined by 

P( Vi = n ) = u-\l -vT 1 ) 71 - 1 , nen. (3.1.10) 



Let Xk(t) — Zk(t — 1) — Zk(t) > be the absolute size of the jump of the kth particle 
from time t — 1 to t, and let q t = 1 — pt- From (6.5) in [19] . in the stationary process 



P{x k {t) = x) 



[1 - up t )q t 1 



x = 



(3.1.11) 



p t (l - up t )q t 1 (u - l) x {uq t ) x x = 1,2,3,... 



We track the motion of particle ^o(^) in a stationary situation. The initial state is 
defined by setting zq{0) = and by letting {^(0)} be i.i.d. with common distribution 



(3.1.10). With k = 0, divide by t in (3.1.9) and take t — > oo. Apply laws of large numbers 



inside the braces in (3.1.9), with some simple estimation to pass the limit through the 



infimum, to find the average speed of the tagged particle: 



lim \z (t) = sup{ux - g(x)} = f(u). 



(3.1.12) 
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For further details please refer to the proof of (7.15) in 

The last equality defines the speed / as / = g + , the monotone conjugate of g. It is 
natural to set f(u) = for u G [0, 1), f(b~ l ) = ), and f(u) = oo for u > b~ l . 

By [16, Thm. 12.4] 



g(x) = sup{xu — g + (u)} = sup {xu — /(«)}• 

«>0 1<w<1/6 



(3.1.13) 



Since Zo(i) is a sum of jumps Xo(k) with distribution (3.1.11) we have the second 



moment bound sup teN E[(t 1 Zo(t)) 2 ) < oo, and consequently the limit in (3.1.12) holds 
also in expectation. From this 



/(«) = - lim E[t~ x zo(t)] = lim Elf 1 x (k)] = E[ar o (0)l 



t— >oo 



k = l 



E^x(u- l) x (uq)- x p(l -up)(l -pY 1 = E 

X=l 



pit (it — 1) 
1 — up 



(3.1.14) 



Next we will find the explicit expression of g(x) from (3.1.14) and (3.1.13). To 
find the supremum of xu — f(x), we compute its first derivative and find it equal to 



x-E 



(b-p) 2 



When E 



< x < E 



(6-p)3 



the equation 



x + 1 = E 



1 — p 



(1 - u pf 



(3.1.15) 



has a solution u G [1, so g(x) = xuq — f(u ). If x < E j^- , then g(a;) = x. If 
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x > E 



i-p 
(b- P y 



/(f). Therefore, 



x/b-b-^l - 6)E 



(6-p) 



c/(a;) = < 



ugE 



p(i-p) 

(l-UQP) 2 



x > b 2 E 
E , 



(i-p) _ i 
(fe-p) 2 1 

< a; < 6 2 E 



(i-p) 

(b-p) 2 



- 1 



< x < E 



i-p 



where «o G (1, 6 1 ) is uniquely defined by the equation 



(3.1.16) 



x+ 1 = E 



(1 -p)(l-u p)' 



Then we need to find the inverse function h(x) = g _1 (x) and then ty^(x,y) = 
yh(x/y). g(x) has three different cases when x takes different values. The first and 
last cases g _1 (:r)can be calculated directly, and for the second case we only need to 
interchange the positions of x and g(x) in their defining equations. Therefore 



h(x) = g 1 (x) 



bx + E 



(i-fr)p 

b—p 



E 



x 



(i-p) 

(1-M0P) 2 



x > E 
E 



p(i-p) 
(fe-p) 2 



i-p 
< x < E 



< x < E 



p(i-p) 

(b-p) 2 



I-P 



with 



r ^oP(i -p) 

:i-%p) 2 



x = E 



(3.1.17) 



Since ^-^(x, y) = yh[xjy\ (3.1.3) proved. 
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To prove (3.1.5) we return to the duality (3.1.13) and write 



g(x) > sup {xu — f(u)} for f(u) = — | —p. 



1<u<1/6 



1-ub 



(3.1.18) 



f'(u) = x is solved by u* = b' 1 (l - ^ . . 
When x > J^r, we have u* G [1, |), and then 



1 , 

(x) > xu* - f(u*) = p(v (1- b)p- \/bx+p) 



Consequently 



g < - {\lb 2 x + a/(1 — b)p) — - = bx — p + 2a/(1 — b)px. 



When x < t£t, the supremum in (3.1.18) is attained at u = 1, and in this case 



g < x < bx + 2a/ (1 — b)px. 



The bound (3.1.5) now follows from ^_^(x,y) = yg 1 {x/y). 



□ 



Proof of (3.1.4) and (3.1.6). The scheme is the same as above, so we omit some more 



details. The inverse of the last-passage time is now defined 



r((a, s), k, t) = min{/ G Z+ : T t ((a, s + 1), (a + I, t + 1)) > k}. 



When k = we define T((a, s), 0, t) = 0. Vertical distance t — s allows for at most t — s 
points with value 1, so the above quantity must be set equal to oo for k > t — s. The 
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particle process {z(t) : t G Z + } is defined by the same formula (3.1.9) as before but it is 
qualitatively different. The particles still jump to the left, but the ordering rule is now 
Zk(t) < Zk+i(t) so particles are allowed to sit on top of each other. Well-definedness of 
the dynamics needs no further restrictions on admissible particle configurations because 



the minimum in (3.1.9) only considers i G {k — t, . . . , k} so it is well-defined for all initial 
configurations {zi(0) : i G Z} such that Zj(0) < Zj + i(0). 

The following can be checked. Under a fixed environment {pj}, the gap process 
{r)i(t) = Zi+i(t) — Zi(t) : i G Z} has i.i.d. geometric invariant distributions P{r]k — n) — 
(i+«)(i+ii) n ' n ^ m dexed by the mean u G M+. In this stationary situation the 
successive jumps x^if) = z k (t — 1) — Zk(t) of a tagged particle have distribution 



l 



P(x k (t) = y) = I 



l+upt 



y = 



From here the analysis proceeds the same way as for the other model. The speed function 
is defined by 



/(«) = - lim E[n _1 z (n)] = E[x (0)] = u(u + 1)E 



P 



. 1 + up 



We now calculate g(x) = sup u>0 {xw — f(u)}. 

When p < x < 1, there exists a non-negative solution uq to the equation 



x = f( x ) = 1 — E 



1 — p 

1 + up)' 



(3.1.19) 
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and it follows that 



g(x) = xu - /(ito) = E 



ugp(l - p) 



L (1 + UqP) 2 

When x < p, we can check the sup is taken at u = and thus g(x) 
g(x) = +oo so 

+oo X > 1 



0. If x > 1, 



^( 2 )P(i-p) 

(1+ M gp)2 



and it has an inverse function 



P < x < 1 
< x < p 



(3.1.20) 



1 — E 



9 [x) 



(i-p) 

(1+mop) 2 



< x < E 
x > 2 i=2 

P 



1-p 

p 



with u defined implicitly in 



x = E 



ulp(l -p) 

(1 + UqP) 2 



Since ty^(x,y) = yg (^), (3.1.4) follows easily. 



To prove (3.1.6), note that 



g(x) = sup{xu — f(u)} > supjx-u — — - — _ \ 

u>0 u>ot 1 + Up J 
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We used Jensen's inequality and concavity of p i-> j^r- From this 



p — px + 2 yp(l — p)x < x < ^-r 2 



and (3.1.6) follows. □ 



3.2 Limiting shapes for exponential models 
3.2.1 Estimate of V G (a, 1) 

In this section, we consider both cases ^(1, a) and ^f(a, 1) for the exponential model 
where some (partially) explicit calculation is possible. 

Let {£j}jez+ be an i.i.d. sequence of random variables < c < £j with common 
distribution m. We assume c is the exact lower bound: m[c,c + e) > for each e > 
0. The distribution function of exponential distribution with random parameter £j is 
Gj(x) = 1 — e~^ x for x > 0. Its mean is h and the variance is 4-. Then the essential 
supremum of the random mean is /i* = c _1 . 

We assume X(i,j) ~ Gj( x ) f° r (^j) ^+ an d write for the limiting time 
constant defined in ( |2.2.4 ). Define fia = f\ coo \ \ m (dQ and a 2 G = Jj coo -j |^m(d£). 



An implicit description of the limit shape was derived in [20J by way of studying an 
exclusion process with random jump rates attached to particles. We recall the result 



from [20] here. One explicit shape is needed for the proof of Theorem 4.1.1 also, so this 
result will serve there too. 

Define first a critical value u* = f< COQ j ]rz^ m (d(,) G (0, oo]. For < u < u* define 
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a = a{u) implicitly by 

u= -^—m(d£). 

J[c,oo) K — a 

The function a(u) is strictly increasing, strictly concave, continuously differentiable and 
one-to-one from < u < u* onto < a < c. We let a{u) = c for u > u* if u* < oo. 
Then define g : M + — > 1R + by 

g(y) = sup{-yu + a(u)}, y>0. (3.2.1) 

The function g is monotone decreasing, continuous, and g{y) = for y > a'(0+) = 
1/ ' hg- It is the level curve of the time constant. The equations connecting the two are 
g(y) = inf {x > : ^g(%, y) > 1} and 

* G (a;, y) = inf {t > : tg{y/t) > x}. (3.2.2) 

We first derive a result of the asymptotic behavior of ^(x,y) close to the y— axis. 
From homogeneity, we can focus on the univariate funtion ^(1, at). 

Theorem 3.2.1. For the random exponential distributions defined above, as a \ 



tf G (a, 1) = Hg + 2a G Va + 0(a). 



Proof. Recall the definition of the limit shape \1/ G (a;, 1) from (3.2.2). From ( |3.2.1 ) one 



can read that tg(l/t) is nondecreasing in t. Thus by (3.2.2) \I/ G (a:, I) = t = t(a) such 
that tg(l/t) = a. 

Next we argue that when a is close enough to 0, <?(!/£) = —Uo/t + a{uo) for some 
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< u < u* with a'(u ) = 1/t. Since a(0) = and a(u*— ) = c, strict concavity gives for 
< u < u* 



c,oo) 



3) (f - «(«)) 



<a'(0+) = {[ C'midO] 

V[c,oo) - 1 



-1 _ 1 



On the other hand, < ^c( a , 1) — A*g < C\/a + Cft where the second inequality 



comes from comparing {Gj} in (4.1.3) with identically zero weights. Thus when a is 
small enough, 1/t is in the range of a'. Consequently there exists uq G (0,u*) such that 
a'(u ) = 1/t, or equivalently, 



c,oo) (f - a K)) 



-m(dt;) = t. 



From the choice of t, a — tg{l/t) = t(—u /t + a(u )) = —u + ta(u ) and so 



(3.2.3) 



tt G (a,l) =t 



ft Uq 
+ 



ft 



+ 



a(w ) a(u ) o(m q ) 7 [Ci0o) £ - a(u ) 



-m(dO- (3-2.4) 



Combining (3.2.3) and (3.2.4) gives 



a = a(uo)" 



coo) (f - «(«o)) 2 



(3.2.5) 



From this 

a(«o) 2 o"G = a(it ) 
Hence we have < a(u ) < yfa/ac- 



1 1 j^m(d£) < ft. 

[c,oo) S 
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When a and hence a(w ) is small, (3.2.5) and the last bound on a(u ) yield 



f I 1 

^ ~ °g = / b — - « W d O 

a(«o) 2 y [c ,oo) (^-o(Mo)) 2 ^ 2J 



= / 2 t ( ; ,o) (3-2-6) 

< 2 / 77 fl( "? } ^ m(d£) = 0(V5)- 



Consequently 



-a 2 



Qi a(«o) 2 G r\l l~\ 
G G = — =— = O(v^). 



Now we put all the above together to prove the lemma. 



*g(«, 1) - //g - 2o- G A/a 

a /" 1 r - 



a(« ) 



a(«o) 7[ c ,oo) £ - a(«o) 
/* r 1 1 

+ / 7 + 7^ a(tto) + O(a(M ) 2 ) - \l g - 2a Gy /a 

J [c,oo) S S 

^ - a G ) + a G a(. ) (a G - j + a ■ O (— j 
= 0(a) as a 4- 0. □ 

We may also approach this problem from a different point of view. It uses the 
following theorem proved in 



Theorem 3.2.2. Two -/M/l queues in series are interchangeable: for any common 
input process A(t), the output processes o/./Mi/l — > /M 2 /l and ./M 2 /l — > /M\/l have 
the same distribution. 



In terms of the language in our last-passage models, the above theorem simply implies 
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that if we pick a realization of the exponential rates the distribution of T(m,n) 

will not be affected if we exchange ^ and £j 2 for < ji < j 2 < n. We can start from 



this idea and reprove Theorem 3.2.1 



Proof. We first look at the case where the distribution m of the exponential rates has a 
finite state space {a\, a 2 , . . . ok}- Assume m({ak}) = Xk, k — 1,2, ... K with J2k x k — 1- 
We now try to approximate ^c{ a i 1) by calculating T( \_Na\ , N) for very large N. 

Fix an arbitrarily small 5 > 0. It is standard result that as iV grows, the number of 
rows that have exponential rates at is in [ [iV(l — S)xk\, \_N(1 + 5)a;fcJ] with probability 
converging to 1 exponentially fast. 

Now we temporarily assume there are \_N(1 + 5)xk\ rows with rate for each k. 
Next we calculate the last-passage time from origin to ( [Na\ , J2k + 5)xu\ — 1). 
Because exponential variables are positive, this result will be no smaller than the actual 
T( \_Na\ , N) if is large enough. 



From Theorem |3.2.2| we can rearrange the rows without changing the distribution of 
the last-passage time. So without loss of generality we let the first [N(l + S)xi\ rows 
have rates a±, the next [N(l + 5)x2\ rows have rates 02, and so on. 

Now we have divided the first quadrant into K horizontal strips, each of which 
has underlying distribution exp(afc). We select the optimal path tt* from the origin 
to ( [Na\ , [N(l + S)xk\ — 1), and record the points at which it exits each strip: 
(M x - I, [N(l + 5)xi\ - 1), (Mt + M 2 - 1, [N(l + 5)xi\ + [N(l + 6)x 2 \ - 1), . . . 
([Na\, EfcLA^ 1 + 5)x k \ - 1). Note M x + M 2 + • • • + M K = [Na\ + 1. 

Now we fix a small e > 0, and define Tk = Then < r&, and ^2 k rk < 
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We look at the intersection of it and the k-th strip. The weight of this portion of 
7r is no more than the last-passage time from the origin to (\_Nr k \, \_N(1 + 5)xk\), with 
weights i.i.d. from exp(afc). Now Theorem 1.6 from [TO] can be applied and it shows 
that the maximal weight in this portion is bounded from above by iV[^(y(I + 5)x k + 
^Jrk) 2 + 5\ with probability exponentially close to 1 in N . In fact for each k the rate 
function 

lim 1 log P (\ ±-T k { [Nr k J ,[N(l + S)x k \)--( ^{l + 5)x k + ^) 2 \ > s) 

depends on the random variable r k . However, {ri, . . . , tk} has a finite state space 
{e, 2e, . . . , rf~|e}, so we still have a deterministic upper bound for the rate function. 

So now we connect all strips and see that with probability converging to 1 exponen- 
tially fast, T( [Na\ , N) is bounded from above by N Y, k [^(v / ( 1 + s ) x k + V^f + <*] • 
By Cauchy-Schwarz inequality, 



^[-(^(l + ^ + v^) 2 

V —(1 + 5)x k + 2 V — + 5)x k r k + V 
k fe k 



'A; 



< 5^ — (i + s)x k + 2 -2^ 1 + 6 ) x * -J2 rk + Yl 



r_k 

a k ' \/ Z T*' a £ ' ^T*' ^ c 



(3.2.7) 



V — (1 + S)x k + 2, / V -^(1 + 5)x fc Va + ATe + -(a + Ke). 



This shows that ^/g(« ; 1) = li m iV->oo ^( |_^ a J j N)/N is P — a.s. bounded from above by 



n, i / /it r 



fc 



-a 
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since we can make 5 and e arbitrarily small. 

Then we can run a similar argument to find the lower bound: first assume there are 
[N(l — 5)x k \ rows with rates a k for each k and rearrange them to get K strips. We 
then pick a path ir e II( [Na\ , J2 k \_N(1 — S)x k \ — 1) in which we let 

Efcl 1 -o)x k a k 

be the number of horizontal movements ix makes in the fc-th strip, and in each strip 7r 
chooses the path that gives the maximal weight. 

With probability converging to 1 exponentially fast in N, the weight of n and hence 



T{ [Na\ , £ fe [N(l - 5)x k \ - 1) is at least N £ fe [±(y/(l - 6)x k + y/M^/N) 2 - 6] . 
From our choice of M k) we get 



iV-s.oo 



V -W{l-5)x k + ^/Ih/Nf 
> V -(1 - + 2 V - v/(l - 5)x k M k /N 

k k 

k ak V * fc 



Then we can use Borel-Cantelli Lemma, let <5 go to 0, and claim the last line above is a 
lower bound of ^g{oi, 1). 

It's worth noting that Yl k ^~ x k is actually the same as /iq, the average of the means 
of the exponential distributions, and J2k ^ Xk should De understood as <7q, the average 



of the variances. Therefore we have already shown that 
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HG + 2a G y/a < # G (a, 1) < pa + 2a G y/a + -a (3.2.8) 

c 



when the distribution m of exponential rates are supported on a finite space. 

We can use discrete distributions to approximate any general distribution m. Let 
{£j,X(i,j)} be a realization of the exponential rates and the weights assigned to all 
lattice points. Choose an arbitrarily small e > 0, couple them with {pj, Y(i, j)} in such 
a way that 

1 

Pj 



and 



Y(i,j) = ^X(i,j). 
Pj 



Then Y(i,j) ~ exp(pj) and p is a random variable with a finite state space. In addition, 

j_ i_ 



we guarantee Y(i,j) < X(i,j) and < ^- — — < e. 



Then we see 



tf G (a,l) = lim -E max X(z) 

n-too n 7ren(|no|,n) z — ' 

> lim — E max y Yu) 

rwoc n 7ren(LnaJ,n) 



> E— +2* E^y/a. 
Po V Po 

Letting 5 go to leads to ^g(®, 1) > Pg + ^gV^ by continuity. 

For the other direction, we can define pj = ^Ji^ £ an d repeat the same argument. 
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This gives the upper bound 



^g(«, 1) < V>G + 2a G Va + -a 



and we have proved the theorem. 



□ 



3.2.2 Estimate of * G (1, a) 



Next, we switch the two coordinates and estimate \I/g(1) ol). Here we see how the tail of 
the random mean /xq creates different orders of magnitude for the a-dependent correction 
term. Qualitative properties of the limit shape depend on the tail of the distribution 
m at c+, and transitions occur where the integrals Jj coo )(£ — c)~ 2 m(<i£) and f< coo \(£ — 
c) _1 m(d£) blow up. ( [20] also addressed this phenomenon.) These same regimes appear 
in our results below. For the case Jj c 00 )(^— c)~ 2 m{d^) = oo we make a precise assumption 
about the tail of the distribution of the random rate: 



The value v = — 1 means that the bottom rate c has probability m{c} = k > 0. Values 
v < — 1 are of course not possible. 

Theorem 3.2.3. For the model with exponential distributions with i.i.d. random rates 
the limit \&g has these asymptotics close to the x-axis. 



3 v G [—1, 1], K > such that lim 




= K. 



(3.2.9) 



Case 1: J\ coo ){£, — c) 2 m(d^) < oo. Then there exists > such that 




for a G [0, «o] ■ 



(3.2.10) 
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Case 2: (3.2.9) holds so that, in particular J, J£ — c)~ 2 m(d£) = oo. Then as 



a \ 0, 



if v E (0, 1] then *g(1, ol) = c 1 + a [ m(d£) + o(a) ; 

J[c,oo) £ ~ c 

if v = then ^g(1, ct) = c -1 — na log a + o(a log a) ; 
ifve [-1,0) then V G (l,a) = (T 1 + Ba^ + o(a^). 



(3.2.11) 

(3.2.12) 
(3.2.13) 



In statement (3.2.13) above B = B(c,k,v) is a constant whose explicit definition 



is in equation (3.2.20) in the proof below. The extreme case v = —1 is the one that 



matches up with Theorem 4.2.3| 



Proof. Equation (3.2.2) gives 



# G (1, a ) = inf{t > : tg(a/t) > 1} = t(a) = t. 



(3.2.14) 



That the infimum is achieved can be seen from (3.2.1) 



Under Case 1 the critical value u* = Jj coo s c(£ — c) 1 m(d£) < oo, and also 



a'{u*-) 



c,oo) 



m{d£) \ > 0. 



By the concavity of a and (3.2.1), for < y < a'(u*—) we have g(y) = —yu* + c. 
Consequently for small enough a 



1 = tg(a/t) = —ac j rn(d^) + ct 

[c,oo) S ^ 
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and equation (3.2.10) follows. 



In Case 2 a'(0+) > a'{u*—) = and hence for small enough a > there exists a 
unique Uo G (0, «*) such that a'(uo) = a/t. Set a = a(ito) G (0, c). As a \ 0, both 
u m* and a /* c. We have the equations 



a'(u ) 1 



£ t 

— —m(d£) = -, 1 = tg(a/t) = -au + ta . 

[coo) (s - fl o) 2 « 



and 



1 , ati 1 
£ = — I = ha 



a AO a ./[c,oo) S — a 



a 



(17, 



c,oo) 



(£ - a o) 



m(df). 



(3.2.15) 



(3.2.16) 



Assuming (3.2.9), start with v G (— 1, 0) U (0, 1). For a small enough e > there are 
constants < K\ < k,2 such that 



Ki(Z-c) v+1 <7n[c,€)<K 2 ((-cy +1 for £ G [c, c + e] (3.2.17) 



and as £ \ we can take /ci, «2 — >■ ft. First we estimate c — do- Fix e > 0. 



a 



c,oo] 



J poo 



2ai 



C+E 



(f - a o 
(£ - «o) 3 



(£ - a ) 3 



for a quantity Ci(e) = 0(s 2 ). The first term above can be bounded above and below 



by (3.2.17), and we develop both bounds together for Kj{ . % — 1 , 2 ? clS 
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rc+e (t _ Y+l 

2^/ \ > d£ + C x (e 



e it - «0) 3 

2Kial I — ^- — " /J d£ + Ci(e 



c+£ [(e-ao)-(c-a ) 1I/+J 



c 

oo 



fc=0 ^ ' 

2K i a 2 A 1/ (c - a o y- 1 - 2 Ki a 2 (" t ^ r^TTT^ " G o)*(c + e - ao)^ 1 + C^e) 

fc=0 ^ 

2K i agA I/ (c-ao)" _1 + Ci(e). 



fe 7 k - v + 1 



(3.2.18) 

Ci(e) changed of course in the last equality. In the next to last equality above we defined 



oo 



v+i\ (-iy 



k J k - v + 1 



Rewrite the above development in the form 

(c - a o y- u = 2nc 2 A u a + a[2A v (KiC% - kc 2 ) + Ci(e)(c - a Q ) x ~ u ]. 

Now choose e = e(a) \ as a \ but slowly enough so that C\(e){c — ao) 1 ^ — > as 
a \ 0. Then also kiOq — > kc 2 and we can write 

c - a = B a^ + o(ai^) (3.2.19) 

with a new constant So = {2k,c 2 A u ) t ^. 
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Now consider the case v G (0, 1) which also guarantees Jj coo \(£ — c) 1 m(c/£) < oo. 



From (3.2.15) and (3.2.19) as a \ 



1 /■ 1 

tf G (l, a) = - + a m(d£) 

a J[c,oo) S — a 



+ a 



1 



C i[c,oo) £ — C 



m(d£) + 0(ai-") + a 



[c,oo) £ — a 



m(d£) — 



[c,oo) £ — c 



+ a; 



^ J [c,oo) S 



m(d£) + 0(a). 



Next the case i/ G (—1, 0). The steps are similar to those above so we can afford to 
be sketchy. 



1 /" 1 

* G (1, «) = — + « / m(d£) 

a J[c,oo) S — a 

1 c-a (c-a ) 2 / "' +; 

= - H n 1 n 1" Q! 



c 2 a 



m[c,0 



d^ + ad(e). 



Again, using (3.2.17) and proceeding as in (3.2.18), we develop an upper and a lower 



bound for the quantity above with distinct constants Ki, i — 1, 2. After bounding m[c, £) 
above and below with c) u+1 in the integral, write (£— c)^ +1 = ((^ — a ) — (c— a )) z/+1 
and expand in power series. 



d£ + aCi(e) 



1 _ 2 1 . 1 

— \- B c a 1 -* + ofa 1 -") + cm* 



OO ✓ 

(c-a o r^( 



i/ + 1\ (-1) ; 



I k — V 



+ a/Ci c-a + £)" V , - — V — +aCi(e) 



Bet 1 -" + o{a 1 - v ) + A Uj2 a(Ki — k)(c — a y + aCi(e) 
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In the last equality the next to last term with the XlfcLo sum was subsumed in the aC\{e) 
term. Then we introduced new constants 



^ = ECT)fc^7 and B = B c- 2 + kB-A^. 

k=0 ^ ' 



(3.2.20) 



As before, by letting e = e(a) \ slowly enough as a \ we can extract \I/g(1, at) = 
c _1 + Ba^v + c^a 1 ^) from the above bounds. 

It remains to treat the cases v = —1,0,1 where integration of the type done in 



(3.2.18) is elementary. We omit the details. □ 
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Chapter 4 



Universality results 



4.1 Limiting shape near the y— axis 



Now we turn to the results on the form of the limit shape at the boundary for a general 
process {Fj}jez+- As explained in the introduction, for ^(a, 1) we find a universal form 



as a \ 0. In addition to the earlier assumptions (2.2.2) and (2.2.3), we need similar 
control of the left tail of the distributions: 



(E[F (x)}) 1/2 dx < oo 



(4.1.1) 



and 



J — c 



esssupF (x) dx < oo. 



(4.1.2) 



Let us point out that ( |2.2.2 ) and (4.1.1) together guarantee E|X(z)| 2 < oo. Let fj,j = 
n{Fj) and cr| = a 2 (Fj) denote the mean and variance of distribution Fj. These are 
random variables under P with expectations \i = E(/i ) and a 2 = E(ctq). Here is our 
main theorem of this chapter: 

Theorem 4.1.1. Assume the process {Fj} is i.i.d. under P, and satisfies tail assump- 
tions d2T2j ) ; fl2~2~3] ) ; pXIfr and gXg . T/ien, as a 4 0, *(a, 1) = /i + 2av^ + o(v^). 
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4.1.1 Proof of Theorem 14.1.1 



We start the proof with the first lemma, which enables us to compare the last-passage 
time limits with different underlying distribution sequences, let {Fj} and {Gj} be er- 
godic sequences of distributions defined on a common probability space under proba- 
bility measure P. In a later step of the proof we need to assume {Fj} i.i.d. Assume 
that both processes {Fj} and {Gj} satisfy the assumptions made in Theorem 



4.1.1 



With some abuse of notation we label the time constants, means, and even random 
weights associated to the processes {Fj} and {Gj} with subscripts F and G. So for 
example fip = E(J xdF (x)). The symbolic subscripts F and G should not be con- 
fused with the random distributions Fj and Gj assigned to the rows of the lattice. 
We write ^ Ber([G(x)-F(x)]+) f° r the limit of a Bernoulli model with weight distributions 
P(X(i,j) = 1) = (Gj(x) — Fj(x)) + = 1 — P(X(i,j) = 0) where x is a fixed parameter. 
An analogous convention will be used for other Bernoulli models along the way. 



Lemma 4. Assume {Fj} and {Gj} satisfy ( p^2[ ) ; (gXTJ) and ( |4L2| ). Then 

for a > 0, 

|tf F (a, 1) - V G (a, 1) - (fi F - fi G )\ 

< 8\/a / (E,\Go(x) — F (x)\) dx + a esssup|F (x) — Gq(x)\ dx. 

J-oo V ' ' J-oo V 

(4.1.3) 



Proof. The right-hand side of (4.1.3) is finite under the assumptions on {Fj} and {Gj}. 
Couple the Fj and Gj distributed weights in a standard way. Let {u{z) : z = E Z+} 
be i.i.d. Uniform (0, 1) random variables. Set Xp(z) = F j " 1 (m(z)), where Fj(u) = 
sup{a; : Fj(x) < u}, and similarly Xa(z) = Gj (u(z)). Write E for expectation over the 
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entire probability space of distributions and weights. 

The following equality will be useful when we relate arbitrary random variables to 
Bernoulli variables: 



X F (z)-X G (z 



I °°{l(X G (z) <x<X F (z)) -l(X F (z) < x < X G (z))\dx. 

J —CO 

Now we compare ^ F (oc, 1) and ^g(oi, 1): 

* F (a, 1) - * G (a, 1) 
= lim — E max V~ " X F (z) — lim — E max V~ " X G (z) 

n-Kx> n 7ren(LanJ,n) n-Hx> n 7reII( [an] ,n) 

< lim — E max S^(X F (z) - X G (z)) 

i r+oo 

= lim -E max V / \ l(X G (z) <x< X F (z)) - l(X F (z) <x< X G (z)) \dx 

ze-rr ^ _oc 

i r+oo . . 

< lim -E / max V<^ l(X G (z) <x< X F (z)) - l(X F (z) <x< X G (z)) \dx. 

We check that Fubini allows us to interchange the integral and the expectation. Since 
F and G are interchangeable it is enough to consider the first indicator function from 
above. Let a be an integer > a. 

/+oo ^ 
-E max V/(l G (z) < x < X F (z))dx 

/+oo i /»+oo 

sup-E max ^ l(X G (z) < x < X F (z)) dx = / *s e r([G(x)-F(x)]+)(a, 1) 

Go(x) - F (x)\ + A^(E\G (x) - F (x)\) 1/2 + aesssup|G (x) - F (x) \ ) dx 



< 



-+oo 

< / (E 



< oo 
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by estimate (3.1.7) and the finiteness of the right-hand side of (4.1.3). Continue from the 



limit above by applying Fubini. Then take the limit inside the ^-integral by dominated 
convergence, justified by the n-uniformity in the bound above. Finally apply again the 



Bernoulli estimate (3.1.7). 



* F (a, 1) - 1 

hoo i 



< lim 



-E max V \l(X G (z) <x< X F (z)) - l(X F (z) <x< X G (z)) } dx 

Tl 7ren(|«n|,n) I J 

00 ZS7T 



/ 1 r 
lim-^E max > l(X G (z) < x < X F (z)) 

+ E max y^(l-l(X F (z) < x < X G (z))) - VlU 

/ + OC 
{^Ber([G( : r)-F(x)] + )(a, 1) + ^ Ber(l-[F(x)-G(x)] + ) («, 1) ~ (1 + «)}c/x 
-oo 

/+oo r 
E(G (x) - F„(z)) + + 1 - E(F (x) - G (x)) + 
oo I 

+ 4v^( ^/e(G (x) - F (x)) + + ^/e(F (x) - G (x)) + ) 

+ a ( ess sup [G (x) — F (x)] + + 1 — ess inf [F (x) — G (x)] + ) — (1 + a) >dx 

\ p F / J 

/+oo /"+00 
a/E|F (x) — Gq(x)| <ix + a / ess sup | Go {%) — Fo(x)\ dx. 
■oo ' J-oo P 

Interchanging F and G in the above inequality gives the bound from the other direction 
and concludes the proof. □ 

For a while we make a convenient assumption that the weights are uniformly bounded, 
so for a constant M < oo, 



P{F (-M) = and F (M) = 1} = 1, 



(4.1.4) 
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then it's easy to see 

a 2 (F ) < M 2 P-a.s (4.1.5) 



and the conditions assumed for Theorem 4.1.1 are trivially satisfied by the uniform 
boundedness. 

Henceforth r = r(a) denotes a positive integer- valued function such that r(a) oo 
as a \ 0. Tile the lattice with 1 x r blocks B r (x,y) = {(x,ry + k) : k = 0, 1, ...,r — 1} 
for (x, y) G Z+. A coarse-grained last-passage model is defined by adding up the weights 
in each block: 

x r ( Z )= X ^)- 

v<=B r (z) 

The distribution of the new weight X r (i,j) on row j G Z + of the rescaled lattice is the 
convolution F r j = F r j * F r j + i * ■ ■ ■ * F r j +r -\. 

We repeat Lemma 4.4 from [13J with a sketch of the argument. 

Lemma 5. Let ^p(x,y) and ^/p r (x,y) be the last-passage time functions obtained by 
using Fj and F r j as the distributions on the jth row, respectively. If r — > oo and 
r-^fa — > as a 10, then 



lim^|^ F (a, 1) - -^ F (ar, 1)1 =0. 
c40 Ja ] r 1 



Proof. Given a path 7r G U(m, nr — 1), consider all the blocks that it intersects; this gives 
a path 7f G Il(m,n — 1) in the rescaled lattice. This path contains almost all the points 
in 7r, with the possible exception at the end point. For example, 7r may contain the point 
(m, nr — 2), but tt does not if r > 2. So there are at most rim + n — 1) — (m + nr — 1) = 
(m — l)(r — 1) points in n but not in 7T, and at most r points in 7r but not in tt. 



(4.1.4) 
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max \ X(z) — max > X r (z) < mrM. 

7ren(m,nr-l) * — ' 7fen(m,n-l) ^ — ' 



Let m = [anrj , divide through by nr, and take limits, finally we arrive at 



lim-^=|*(a,l) - -^(ar,l)\ <YmvMr^a = (4.1.6) 

aiO Jet V c40 



□ 

Let fi r)V and V r>y be the mean and variance of F ry : 

r-l r-l 

^r,y = ^ry+i, and = °"ry+r 

i=0 i=0 

Let <& r)3/ be the distribution function of the normal A/"(/i r>2/ , distribution, and $ ri?/ 
the distribution function of J\f(rfiF,V rty ). The difference between $ r ^ y and is that 
the latter has a non-random mean. We shall also find it convenient to use {Xj} as a 
sequence of independent variables with (random) distributions Xj ~ Fj. For the next 
lemma we need to assume {Fj} an i.i.d. sequence under P. 

As in [13], a key step in the proof is the replacement of the rescaled weights with 
Gaussian weights, which is undertaken in the next lemma. 

Lemma 6. Assume {Fj} i.i.d. under P. If r — > oo and rJ~a — > as a 1 ; then 



lim J— |tf Fr (ar, 1) - *$ r (ar, 1)| = 0. (4.1.7) 



c40 TV a 



Remark 7. TTie following proof contains many inequalities where we need to use letters 
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to denote proper constants. For simplicity all constants that does not depend on r, a and 
M are subsumed in a single notation C . 



Technically we could also subsume M into C and the proof of Theorem J^/LA is not 



affected. However, we write M explicitly to help us prove the next theorem 4-1-2 



Proof. We will use (4.1.3) for the processes {F rj y} ye % + and {<& r ,y}yez+ and with a re- 
placed by or. On the right-hand side there are two terms. We will handle the second 
term first. 

To estimate the second integral, we discuss over the value of x. For x > 2rM, we note 



that F riV (x) = 1 since we assume (4.1.4). Now we turn to & r ^ y (x). We quote Theorem 
1.4 from [1] and get 



$ r Jx)=l- , 1 exp(- (s tlr ' y)2 )ds 



2 



>1 _ V 7 ^ exp (- (x ~^ ) ) 

V27r(x - /I rt y) 2V r: y 



VrM 2 . (x-rM) 2 , 

>1 == exp — 

v / 2^(2rM -rM) K 2rM 2 ' 



(4.1.8) 



This gives 



1 , (x-rM) 2 . 



\F r Jx) - <S>,, v (x)\ < -jL= exp(- (X 2r ^ f)2 ) (4.1.9) 



for x > 2rM. From symmetry a similar inequality 



\F r , y (x) -* r *(*)\ < 7^ eX P(- (X 2 tl^ ) ) (4L10) 



holds for x < —2rM. 
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Now we are ready to claim 



-oo 



ar / esssup|F ri0 (x) — <& r fl(x)\ dx 



oo ir 

^ , /- 2rM 1 / (x + rMf ,, M f 2rM 1 . 
<ar< / exp — — — ax + / 1 • dx 

/•+ 00 1 , (s-rM)\ , 

+ 1m 71S eXP( -^^ )i 



<Car{ M + 4rM ^ < CMar 



,2 



Next we estimate the first term on the right hand side of (4.1.3 ). We will need Theo- 
rem 5.17 of [14], which states that if independent mean random variables Xi, X 2 , X 3 , . . . 
all have finite third moments, then they satisfy the estimate 



i=l 



where B r = X^i=i Var(Jfj), $ is the standard normal distribution function, and A is a 
constant that is independent of the distribution functions of X%, X%, . . . , X r . 

Recall that we assume {Fj} i.i.d. under P, so <7 2 (Fj) are i.i.d. random variables. For 
an arbitrary < e < 1, say e = |, we define C/ r as the event £)£~ ofu+i — r (l ~~ £)Eo"q = 
IrEcTg. Then it is standard result that F(U r ) converges to 1 exponentially fast as r goes 
to infinity. 

With probability F(U r ) we get 

iF r , tf (x) - * r M\ < ^ E:: ^ r ^~7r |3 (i + v-y*\x - ^, y \r 3 

\l^i=0 a ry+i) (4.1.11) 

< + M - 1 r- 1 ^\x-iir,y\)~ a 
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where the second inequality used the assumptions -P(|Xj| < M) = 1 and the property 
of U r , and C is a proper constant that depends on A and Eo-q. 

Next we note this trick using Cauchy-Schwarz inequality: for a probability density 
/ on R and a function H > 0, 



For the calculation below take 5 > and f(x) = Ci(l + |rr — r / u F | 1+<5 )~ 1 for the 
right constant c\ = Ci(5) to make f(x)dx = 1. Again factors that depend on 5 are 
subsumed in a constant C in each of the following steps. 

Over the event U r , 

r+oo -y 




1/2 



Then we get 




(4.1.12) 




by a change of variables x = fi r fi + yM\fr 



(4.1.13) 




< CM 4 (|/i ri0 - rfM F \ 1+5 + m 1+ V 1+ ^ 2 ). 



Over the event U£ we use (4.1.9) and (4.1.10) to bound the integral 
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-+oo y 

\F r o(x) - $ r0 (x)\dx 

i i r I 

' — oo 



< 



/(*) 



We use a change of variables y = and the first term above 

r-Vr m v 2 

= (l + \y/rMy -rM -r/j F \ 1+s )^=exp(-^-)dy 

J-oo v27r 2 

/•"V? /If 7 .2 

<C / (| v / ^M?/| 1+ ' 5 + |rM + r/i F | 1+<5 )^exp(-^-)(i|/ 
J-oo v27T 2 

<C f ^M 1+ V^M 1+5 ^Lexp(-^)^ 

./-oo v27T 2 

+ C / r i+*M 1+ *^Lexp(-^)dy 
./-oo v27r 2 

< C ( r (l+5)/2 M 2+5 + r l+<5 M 2+<5) < Cr 1+8 M 2+S 

The third term follow the same upper bounds. 

The second term is simply bounded by C ■ rM ■ (rM) 1+5 = CM 2+S r 2+S . Then 



/+oo 1 
— |F r , (x) - ® r ,o(x)\dx < CM 2+S r 2+5 . (4.1.14) 
-oo J\ X ) 



Continue from (4.1.12), and keep in mind that P(U£) decays exponentially as r grows, 
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we have 



. r + °° / \V2 

far I ( E\F rfi (x) - $ r ,o(^)| ) dx 

J —oo 

[+co 1 

< V«H E / — - 1 F ri0 (x) - $r,o (a;) I da; 

1 i-oo /W 



+' x l s 1/2 

(4.1.15) 

1/2 



< CV^{M 4 E|^ n0 - r/i F | 1+5 + M 5+ V 1+5)/2 + M 2 +V +5 P(f/ r c )} 

< CM (5+5)/2 a 1/2 r (3+<5)/4 . 

Here we used the fact that /i r — jj, F is a sum of independent bounded mean-zero 
variables, so 

/ \ ( 1 +< 5 )/ 2 

E|Ai r , - r/i F | 1+ < \E\n r! o - v\i F | 2 J 

= [rE(// - /xf) 2 ] {1+5)/2 = CM 1+5 A 1+5 " 2 . 



To summarize, with these estimates and (4.1.3) we have 

' |*,. r (ar,l)-** r (ar,l)| 



<_^( M (5+<5)/2 a l/2 r (3+5)/4 + Mm>2) ( 4L16 ) 



=C(M (5+^)/2 r (-l+5)/4 + Mrv ^ ) _ 

By choosing 5 < 1, then assumptions r — > oo and r-y/a — > make this vanish as a — >■ 0. 
The proof is completed. □ 

Remark 8. The proof presented above distinguished the two events U r and U£ because 
(4.1.11) only works when V rjV = YliZo a j can ^ e bounded away from zero. Therefore on 



US: we used a different approach. 
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We can actually find an alternative proof of this lemma. In addition to (4.1.4) let's 



also assume that variances are uniformly bounded away from zero, i. e. for a constant 
< Co < oo, 

P{<t 2 (F ) > c } = 1. (4.1.17) 



Note that then cq < a (Fq) < M . A direct consequence of this is that now (4.1.11) and 



(4.1.13) holds F-a.s., so 



/+°° / \ 1/2 

(E\F rfi (x) - $r,o( x )\) dx < Ca 1/2 r (3+5)/4 . 
oo 



Note that in equations above and below we subsume M into C for simplicity since we 



will not come back to them any more. For the second term on the right in (4.1.3) 7 



nr I esssup|F ri0 (s) — ® r ,o( x ) \ dx < Car 1 / 2 j ess sup ( 1 + 

-oo P 

<Car 1 ' 2 ' 



\x — fJ,, 



r,y\ 



-3 



dx 



- ,U/ -rM -x\~ 3 , f rM , f +co f x-rM\~ 3 , 
1 H — = ] dx + / cfe + / ( 1 + — = dx 



My/r J J_ rM J rM V J 

{POO poo 
My/r J u~ 3 du + 2rM + M^fr J u' 3 du 

Car l/2 (M^/r + 2rM) < Car 3/2 . 



Therefore we see (4.1.16) still holds and Lemma^is proved under (4.1.17). However, 
we eventually we have to show Theorem 4-1-1 without ( |4.1.17 ), so then we try to lift this 
assumption. 

For e > 0, let {W(z)} be i.i.d weights with distribution H defined by P(W(z) = 
±e) = 1/2. Let Fj = Fj*H be the distribution of the weight X(i, j) = X(i, j) +W(i, j) . 
Let and be the time constants of the last-passage models with weights {W(z)} and 
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{X(z)}, respectively. The Bernoulli bound (3.1.7) gives the estimate ^h{ x iV) < 
The corresponding last-passage times satisfy 



Tdz) - T H (z) < T F (z) < Tdz) + T H (z 



where Th{z) uses the weights —W(z). In the limit 



*p(a, 1) - ^e^fa < V F (ae, 1) < *s(a, 1) + 4ey/a. 



(4.1.18) 



Since a 2 (Fj) = o- 2 (Fj) + e 2 , {Fj} satisfies ( |4.1.17p . Once {Fj} j&+ satisfies fl4.1.25| ), 
t/ien so does {ij-}jez + > because fip = \ip and e > can 6e arbitrarily small. 



After the discussion of Lemma [6] we make a further approximation that puts us in 
the situation where all sites have normal variables with the same mean. 



Lemma 9. Let \E , $ r . and be defined as before, and again r^fa — > as a — > 0. Then 



lim — =|¥* r (ar, 1) - (or, 1)1=0. 



Proo/. For z = (ij) e Z 2 + , let X {t \z) have distribution $ rj so that X^(z) = X^(z) - 
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fj, r j + rfj,p has distribution $ r j. Now estimate: 



% (ar.l) = lim - max y^X {r) (z) 

* r n^oo n vrGn(LanrJ,n) 

ZG7T 

< lim — max X^ r \z) + lim — max Y^f— /V 7 + t^if) 

n-*oo n 7rGn(LonrJ,n) n-5-oo n 7reII( [anrj ,n) 

1 n 1 

< ^ (ar, 1) + lim — } (—Mr 7 + t ^f) + lim —2Mr ■ \anr\ 

n-s-oo n v ' n-5>oo n 

i=o 

= $$>r, 1) + IMar 2 . 

Note that in the second to last step we used the fact that when is assigned with 
—fXrj+rfXp, every admissible path in n[|_a:nrj , n) contains at least one of each —{i r j+rfi F 
for every j — 0, 1, . . . , n. Their average converges to by the law of large numbers. There 
are also [anr\ additional points all bounded from above by 2M, which contribute to 
the third term. 

The opposite bound (ar, 1) > ^^ r (ar, 1) — IMar 2 comes similarly. So 



(ar,l) -** r (or,l)| < 2Mr^/a~ (4.1.19) 



r^ 1 *' 

and the lemma follows. □ 

Let us separate the mean by letting ^ T)V denote the N(0, YaZo a %+i) distribution 
function. Since the last-passage functions of the normal distributions satisfy 

ttx (ar, 1) = r/i F (l + ar) + %,)(w, 1), (4.1.20) 



we can summarize the effect of the last three lemmas as follows. 
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Lemma 10. Assume {Fj} i.i.d. under P, and assume r = r(a) satisfies r — > oo and 



r^fa — y as a \. 0. Under assumptions (4.1.4) 



]im^=\V F {a, 1)-Hf- -* s w(ar, 1)| = 0. 



c40 ^/a 



(4.1.21) 



In order to deduce a limit from (4.1.21) we utilize the explicitly computable case of 



exponential distributions from [20J, and use the results proved in Chapter 3.2 We need to 



match up the random variances of the exponentials with the variances crj of the sequence 
{Fj}. Thus, given the i.i.d. sequence of quenched variances cr| = o- 2 (Fj) that we have 



worked with up to now under condition (4.1.5), let £j = 1/cr,- and Gj(x) = 1 — e ^ x the 
rate £ 3 - exponential distribution. Then is an i-i-d- sequence of random variables 



£j > with distribution m. Since we assume (4.1.4), the sequence {£j}jez + is bounded 
away from zero. We can assume c is the exact lower bound: m[c, c + e) > for each 
e > 0. Gj has mean and variance fJ>(Gj) = l;^ 1 and a 2 {Gj) = ^j 2 = a 2 . 



Assumptions (2.2.2) and (2.2.3) are easily checked, and so the last-passage function 



\&g is well-defined. We would like to apply Lemma 10 to this exponential model, but 



obviously assumption (4.1.4) is not satisfied. To get around this difficulty we do the 



following approximation which leaves the quenched means and variances intact. We 
learned this trick from 1 1 3 1 . 



Let Yj denote a G^-distributed random variable. For a fixed r > 0, let 



E(Xj\Yj > r) and Wj = E(Y 2 \Yj > r). 
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The quantities 



(rrij — r) 1 



(rrij — r) 2 + Wj — m 2 



and Uj = — 



T 



satisfy the equations 



[1 — Sj)r + SjUj = rrij and (1 — Sj)r + SjUa = Wj. 



Then < Sj < 1, Uj > r and Wj > t 2 . Define distribution functions 



Gj(x) = ( 



Gj(x) < x < t 

1 — Sj[l — Gjij)} T < X < Uj 



1 



X > Mi 



(4.1.22) 



Yj ~ Gj satisfies EYj = EYj and EY 2 = EY 2 . Moreover, for any fixed r > 



u ; 



E(Y 2 \Y 3 > t) - t 2 



2 2 

T = + T < 1- T 

^•1^ >r)-r Pi ~ c ' 



so {Gj} are all supported on the non-random bounded interval [0, 2/c+r]. Consequently 



Lemma 10 applies to G. We can draw the same conclusion for G once we have the next 
estimate: 



Lemma 11. Given e > 0, we can select t large enough and define Gj as in (4.1.22) so 
that 

Iim^|^ G (a, 1) - ^(a, 1)1 < e. 
c40 J a 
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Proof. This comes from an application of Lemma 4 Gj = Gj on (— oo, r) and 1 — Gj < 



1 - G 3 on all of E. 



GV"' 

"+°° / ~ \ 1/2 /•+°° _ 

< 8-v/a / [ M\Gq(x) — Go(x)\ ) dx + a ess sup | Go (x) — Go(x)\ dx 

-oo J J~oo P 

"+°° / : \ 1/2 

< 8Va / (E\1-G {x)\+E\1-G {x)\) dx 



+ a (esssup|l — Go(x)\ + esssup|l — Go(x)\)dx 

r+°° , xi/2 r+°° 

<8V2aJ (E\l-G (x)\j dx + 2aj esssup|l - G (x)\dx 

/ +OQ cx f +OQ 

exp(— — )dx + 2a J exp(— cx) dx 



16y2a: ct 2a , . 

ex P(~v) + ~~ exp(-cr) 

c 2 c 



(4.1.23) 



Now we see lim a ; ^\^g(oc, 1) — \E , g(a, 1)| can be made arbitrarily small by choosing 
r large. □ 



So Lemma 10 and Lemma 11 together show that 



lim -Utt G (a, 1) - Ea - -% r (ar, 1)| = 0. 
<x|.o A/a r 



(4.1.24) 



It remains to perform an explicit calculation on ^c^a,!). As before, utilize the 



notation \iq = E£ 1 and a 2 G = E£ 2 . In Theorem 3.2.1 We have already computed that 



lim -^=\^ G (a, 1) - no - 2o G ^[a\ = 0. 
c40 wa 
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This result combined with (4.1.24) gives 



lim — = I -^5 (or, 1) - 2a G y/a\ 
c40 y/a 1 r 1 



Substitute this back into (4.1.21) and recall that a F = o G . The conclusion we get is 



lim —=\^p{oi, 1) — llf — IofxToA = 0. 
c40 Ja 



(4.1.25) 



So far we have proved Theorem 4.1.1 under the assumption (4.1.4). As the last item 



of the proof of Theorem 4.1.1 we remove this uniform boundedness assumption. Suppose 



{Fj} satisfy the conditions required for Theorem 4.1.1, but there is no common bounded 



support. For a fixed M > define the truncated distributions 



x > M 



Fj, M (x) = { Fj ( x ) _ M < X <M 



x < -M. 



Let [Am, &m an d ^F M (%,y) be quantities associated to {Fj >M }. 



From (4.1.3) and the conditions assumed in Theorem 4.1.1 



r+°° / n 1/2 r+°° 

<8 / (E|F (x) — Fo t M(x)\ ) dx + y/a / esssup|F (o:) — F 0j m(%)\ dx 

J-oo V ' J-oo P 

- M / \l/2 /- 00 / Nl/2 ■ 

(E|F (x)|) dz+ / (E|l-F (x)|) 

M r+oo 

esssup|F (x)| dx + / esssup|l — F (x)| cix < e. 

D P JM P 



(4.1.26) 



+ v a 
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The last inequality comes from choosing M large enough, and is valid for all a < 1. 
Since E(£'X 2 (0,0)) < oo, dominated convergence gives om — > o~ and so we can pick M 
so that \a — <tm\ < £• Now 



'^fpfa, 1) — a — 2ay/a\ < —=\^fp M (a, 1) — llm — 2&mV®\ + 2e. 
a \ a 



Since e is arbitrary and limit (4.1.25) holds for {-Fj,m}, we get the conclusion for the 



sequence {Fj}. This concludes the proof of Theorem 4.1.1 



4.1.2 An improvement on the error o{^/a) 



The error term o(yfa) can still be improved. As was shown in Theorem 3.2.1, in the 



exponential model 0(a) is a more accurate estimate. In this section we will take a closer 



look at the approximations we used along the proof of Theorem 4.1.1 and analyze the 
order of a. 



Theorem 4.1.2. Assume the process {Fj} satisfies the assumptions in Theorem 4. 1.1 



i.e. itisi.i.d. underF, and satisfies tail assumptions (2.2.2), (2.2.3), (4.1.1) and (4.1.2). 



In addition, assume the uniform bound (4.1.4), then for any e > 0, as a \ 



1) = fip + 2opva + o(a 



(4.1.27) 



Proof. We give a list of all approximations in the previous subsection: when {Fj}j £ z. 
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satisfies (14.1.41) 



1 



1 



COb :-=|tf(or, 1) - -tf (or, 1)1 < Mrx/a, 
a r 



(Olfib |* Fr (« r ; !) - 1)1 < C(M (5+<5)/2 r ( - 1+5)/4 + Mr^a) 



1 



(4.1.19) :^^|^$ (or, 1) - ^$ r (ar, 1)| < 2Mr^, and 
1 r rJa 



(4.1.20) (or, 1) = r/x F (l + or) + % r) (ar, 1). 



To sum up, (4.1.21) can be rewritten as 



\V F (a,l) -ti F --%(r)(ar,l)\ < C{M^ +5),2 r^ l+5)/ ^ + Mr a) 



(4.1.28) 



for a proper constant C. 

Now we recall the approximation for the exponential model. For exponential distri- 



butions {Gj}j<=.z + , if we define {G f j}jez+ uniformly bounded by M as in (4.1.22), then 



(4.1.23) gives: 



* G (a, 1) - !)l ^ Cexp(-cM)V^ 



for a proper constant C . 



The above equation, together with (4.1.28) applied to {Gj)jgz + , implies 



|* G (a, 1) - /i G - -%r)(ar, 1)| < C(exp(-cM)ya] + M (5+5)/2 r ( - 1+<5)/4 Va + Mm). 

(4.1.29) 

Here C does not depend on M, r, and a. 
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In the above equation the left-hand side is independent of M, so we make both M 
and r functions of a. Let r(a) = and M(a) = ——In a. Make S small enough, 



then the right hand side of (4.1.29) is bounded by C(— — In a) 5+s a 5 ~ s 
a\0. 



as 



Recall Theorem 3.2.1 we get 



-(r) 



'ar, 1) = 2ac\/a + o(«s" 



(4.1.30) 



Now let us come back to ^(a, 1). If {Fj}j & z + are uniformly bounded, then M is a 



fixed constant and the right hand side of (4.1.28) is just o(as e ), and this shows that 
Ve > 0, 

3 

ty(a, 1) = fiF + 2a F y/a + o(a^ £ ). 

□ 

Remark 12. TTie proof did not treat the case when {Fj}j £ i + is not uniformly bounded. 



The difficulty is that (4.1.26) does not give precise computability on how the right-hand 



side would change according to M . We surely need additional assumptions in this case. 



4.2 Estimates for limiting shape near x-axis 



We turn to the case ^(1, a). As we have seen in 3.2.3, the results will be qualitatively 



different from Theorem |4.1.1| The leading term will be the essential supremum of the 
mean instead of the averaged mean and we will see different orders for the first a- 
dependent correction term. Universality results are desired but much more difficult to 
achieve. In this section we will see some estimates of ^(T, a) under various conditions. 



72 

We first present a result which gives an upper bound of a) in a very general set- 
ting. We will use F as superscripts or subscripts when we want to stress the dependence 
(of a probability, or expectation, etc.) on a distribution function F that is in the state 
space of P. Again, /i(F) is the mean of the distribution F and cr 2 (F) is the variance. 

Theorem 4.2.1. Assume {Fj}j £Z+ is i.i.d. and there exists constants t > and K < 
oo such that 

esssuppSV !*-^! <K . 

Let [/,* = ess sup P [J>(F), then 

^(l,a)=fi* + 0(Ja\og-). (4.2.1) 
V a 

Proof. We will use u> to denote a realization of {Fj}j €li+ , and use P w and E w as the 
corresponding quenched probability and expectation. Recall the definition 

Wl,a)= lim - max X(z). 

n-too n 7ren(n,LnoJ) ^— ' 

For any specific path 7r, we write T(7r) = XLe?!-^- 2 )- We also use as the leftmost 
site traveled by it in the j-th row. Then we derive the following estimate: for any positive 
u and t, 
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P(T(tt) > nu) = EP w (T(vr) > nu) < E(e" ntu E a; (e* TW )) 
= e - n ' u E[exp(^log£ w (e' x(z) ))] 



2£7r 

\na\ Tij+i 

e -ntu J-J Eexp [^ log^(e* x(i ' j) 

\na\ 

E 

j=0 



(4.2.2) 



exp j-ntw + ^ logE[exp((l + m i+1 - rrij) \ogE Fj e tx 



We need the following inequality in order to proceed: if we have a random variable 
Y and positive real numbers t±, ...,t m with J2j tj = ^ then 

5>g£^ y < 5>g[£(e ty )]T = ]T ^log[£(e ty )] = log[£(e ty )]. (4.2.3) 



Continuing from (4.2.2), as n /*■ oo 



P(T(tt) >nu) < exp{-ntu + logEexp(n(l + a) log£ F e tx )} 
= exp{-w[ta- (1 + a) ^_ logEexp(n(l + «) \ogE F e tx )]} (4.2.4) 

<exp|-n[tw - (1 + a)A(t)]}, 

where A(£) = ess sup P A F (t) and A F (t) = \ogE F e tx . 

From the assumption that there is to > and Kq < oo such that E F e to ^ x ~^ F ^ < Kg, 



it follows that for t G (0, to), 
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A F (t) =n{F)t + log(l + ^t 2 + E F Y, -{X- ^F)) k 

k=3 



<»(F)t 



a\F) 



t 2 + K t 



(4.2.5) 



and therefore if we denote a* = esssup P <r(F), 



a 



*2 



A(t) < fft + t 2 + K t 3 , 



from which we obtain 



*2 

»(T(tt) > m) <exp|-n[tM- (1 + a)(fft + °—t 2 + K t 3 )]} 

*2 

= exp{-n[(u- (l + a)//*)t- (l + a)^-t 2 - (1 + a)i^ 3 ]}- 



(4.2.6) 



We take u = fi*(l + a) + e and i = .*2n +a ) > where £ = 4<7*-v/alog -. When a is 



small enough, e is small and we get t G (0, to). Then (4.2.6) becomes 



P(T(7r) > mi) < exp 



-n 



(1- 



2if £ 



2cx* 2 (l + a) v a* 4 (l + a 



< 



exp 



-n 



4er* 2 (l + a) 



(4.2.7) 



when a and thus e is small enough. 

The last-passage time T(n, \na\ ) = max 7rg n( n ,Lnaj) T(n). Here the maximum is taken 
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over a pool of ( n+ [ i nQj ) paths, so we can use Stirling's formula to get 



(n + I na I 
P(T(n, |ria|) >n«) <( ^ J ] oxp 



— n 



4a* 2 (l + a) 



, c 

<— exp 



-n 



C 

(a log a - (1 + a) log(l + a) + 4(T * 2(1 + a ) ) 



(4.2.8) 



Plug in the expression for e we have 



a logo; — (1 + a) log(l + a) + 



4^(1 + a) = TT^ bg a " (1 + a) l0g(1 + ^ 



As a \ 0, the first term on the right-hand side above has order alog^, whereas the 
second term has order a. So when a is small, 



a log a — (1 + a) log(l + a) + 



4a* 2 (l +a) 



> 



and hence 



^^P(T(n, L^«J) > nil) < oo, 
which by Borel-Cantelli lemma gives 



< u = (SCI +a) +4<r*Walog-, 

a 



(4.2.9) 



and we can claim as a \ 0, 



^(l,a)=/i* + 0(Xlog-). 

V ct 



(4.2.10) 
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□ 

The order <J a log ^ should be a rather conservative estimate. In Theorem 3.2.3 and 
in some other settings we find that the first a-dependent term in \I/(l,a) is no more 
than \fa. There has not been a very general statement about the necessary condition 
for the order \fa so far, but next we will see two sufficient conditions. The first one is 
uniform boundedness. 

Theorem 4.2.2. Let {Fj}jez + be an ergodic sequence of distribution functions satisfying 



the conditions listed in Proposition \2.1\ Assume the existence of M > such that 

P{F (-M) = 0, F (M) = 1} = 1. (4.2.11) 

Again define fi* = esssup P n(F), then as a \ 0, 

*(l,a) = n* + 0(y/a). (4.2.12) 

Proof. We first prove an upper bound for ^(1, a). We start by increasing all the weights 
X{z) by moving their means to /i*, so that \l/(l,a) for the shifted weights gets no 
smaller. Then we subtract the common mean fi* from the weights. Therefore we can 
assume fi(F) = for all F. The weights X(z) are still uniformly bounded, and without 



loss of generality we still assume (4.2.11) for the shifted weights with the bounds still 
denoted by M. 

Fix a realization of {Fj}j & % + , and the lattice point zq = (0,0). Let iV be a positive 
integer. According to whether the path 7r goes through z or not, and in case it goes we 
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also separate the weight at z , we write 

max S^X(z)=Av(B + X(zo)) = B + (A-B)VX ( zo), (4.2.13) 

7ren(n,LATaJ) ^— ' 

where A = max Z0 ^ n J2 z eiv X ( z ) and B = max ^enJ2 z eA{M X ^' Botn A and 5 look 
complicated but we only need to treat them as some random variables. Let G(y) denote 
the distribution of A — B, then the quenched expectation 

E[(A-B)VX(zo)]= [ xVydF (x)dG(y) = 1(1 xV y dF {x))dG(y). (4.2.14) 

We now take a closer look at J M M x V ydF (x). The only nontrivial case is when 
y G [— M, M], integration by parts gives 

/M pM 
xVydF (x) = yF (y) + x dF (x) 
■M Jy 

pM 

= yF (y) + [M — yF (y)) - / F (x) dx (4.2.15) 

pM 

I F (x)dx. 

J V 



-M 

= M 

' u 



Next we try to minimize the integral J^ 1 F(x)dx when F(x) is selected from mean 
zero distribution functions supported on [— M, M]. Suppose the value F(y) is given and 
F(y) > |) then obviously J^ 1 F(x)dx > (M — y)F(y), in which the equal sign can be 
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achieved for 







x < -M 



F(x) = { 



M+y 



(^IMM+vm - M <x<y 
y < x < M 
x > M. 



(4.2.16) 



F{V) 



If F(y) is known and F{y) < \ , the above function F{x) would not work since it is then 
not non-decreasing. Now we have 



rM I'M py 

/ F(x)dx = / F(x)dx - i F(x)dx 

Jy J~M J-M 

/M r-y 
xdF(x)- / F(x)d. 
-M J-M 

/y 
F(x)dx 
■M 

>M-F(y)(y + M). 



The equality holds when 







x < -M 



F(x) = < 



F(y) -M <x<y 

(l-F(y))M-yF(y) < M 

M-y il — X ^ 1V1 

1 x> M. 



(4.2.17) 



(4.2.18) 



We now summarize the above two cases and see J M F(x)dx > \{M— y), with equality 



when F(y) = h, which corresponds to the distribution function 
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F M (x) 



x < -M 



I -M < x < M 



1 M < x. 



(4.2.19) 



Notice that F (x) puts half probability on M and — M each and does not depend 



on the value y. Back to (4.2.14), we see that for any random variables A and B, 



E[(A — B) V X(zo)] is maximized as long as we let X[zq) follow F M (x). Running this 
argument for all z G {0, ...,N} x {0, |_iVaJ}, we obtain 



E max > X(z) < E max > X f m(z) 



7ren(n,LVoJ) 



2G7T 



7ren(n,LVaJ) 



ZdlT 



Taking limits and using (1.0.3) gives 



(4.2.20) 



) = 2M v / a + o( x /a). 



(4.2.21) 



If we consider the effect of /i*, we get (4.2.12). 



□ 



The next two theorems also prove that a) is a constant plus order 0(y/a). They 
relax the assumption of uniform boundedness and use the ones from [13J. The fmiteness 
of the state space of P plays an important role in the proofs of both theorems, but it 
does not seem to be a necessary condition for the results. It would be great if we can 
design a different approach and remove this fmiteness condition. 
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Theorem 4.2.3. Assume the process {Fj} of probability distributions is stationary, 
ergodic, and has a state space of finitely many distributions H%, . . . ,Hl each of which 
satisfies Martin's FTStf hypothesis 

poo rO 

/ (I - H^x)) 1 ' 2 dx + / H^x) 1 ' 2 dx < oo. (4.2.22) 

Let fi* = max£fi(He) be the maximal mean of the Hi's. Then there exist constants 
< C\ < c 2 < oo such that, as a I 0, 

fx* + ci V ^ + o(v / ^) < < H* + c 2V / « + o(v / «)- (4.2.23) 



Proof. The lower bound in (4.2.23) can be proved by applying Martin's result (1.0.3) 
to the homogeneous problem where a maximal path is constructed by using only those 
rows j where Fj = H^, the distribution with the maximal mean ii* = fi{Hi*). This is 
fairly straightforward. 



To prove the upper bound in (4.2.23), we again start by increasing all the weights 
X(z) by moving their means to if. Then we subtract the common mean /i* from the 
weights, so that for the proof we can assume that all distributions Hi, ... , H^ have mean 
zero. 

Create the following coupling. Independently of the process {Fj}, let {Xf (z) : 1 < 
£ < L, z G Zi} be a collection of independent weights such that Xi{z) has distribution 
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Hi. Then define the weights used for computing a) by 



L 



Begin with this elementary bound: 



ty(l,a) = lim -E max 

n-+oo n L 7ren(ra,LcmJ) 




(4.2.24) 



< > lim — E max 

^n-Joon L 7ren(n,LanJ) 



^JQ(z)/ {F . =if , } ]. 



The next lemma contains a convexity argument that will remove the indicators from the 
last-passage values above. 

Lemma 13. Let V be a sub-a-field on a probability space (Q, J 7 , P), D an event in V, 
and £ and rj two integrable random variables. Assume that Erj = 0, rj is independent of 
V, and £ and rj are independent conditionally on V. Then E[£ V (tjId) ] < E[£ V rj]. 

Proof. By Jensen's inequality, for any fixed i6K, 



Integrate this against the conditional distribution P(£ e dx \ V) of £, given V, and use 



x V E(r)\V) < E(xVr]\V). 



Since r/ is independent of V and mean zero, 



iV0< E(x Vr)\V). 



the conditional independence of £ and rj: 
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E(£ V0|Z>) < E(£Vri\V). 



Next integrate this over the event D c : 



E[l D c- £V(tiI d )] =E[l D c- £V0] <E[I DC - £Vt?]. 



The corresponding integral over the event D needs no argument. 



□ 



Fix a lattice point Zq = (io,jo) f° r the moment. We split the maximum in (4.2.24) 



like the way we did in (4.2.13): 



max £ X e (z)I {Fj=He} = B+(A-B)V (X e (z )I {Fjo=Hi} ) 

7rGll(n, \na\ ) ' * 



where 



A = max }] Xj(z)I{F i= H t } and B = max 2J X e (z)I {Fj=H 



I' 



2S7r\{2 } 



Now apply Lemma 13 with £ = A — B, r] = X^zq), and -D = {Fj = Hp}. Given F, 



JO ' 



A — B does not look at X^(^o), so the independence assumed in Lemma 13 is satisfied. 
The outcome from that lemma is the inequality 



E 



max Y, X ^ z )hF J =H l} <E[AV(B + X e (z ))} 



This is tantamount to replacing the weight Xe(z )I{ F . =He y at z with Xe(z 



We can repeat this at all lattice points z in (4.2.24). In the end we have an upper 
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bound in terms of homogeneous last-passage values, to which we can apply Martin's 



result (1.0.3): 



#(l,a) < y lim -e[ max V XAz) 1 = V a) 

V ; "^»^oon L 7ren(n,lanJ) ^ V '\ ^-^ £V 

t=\ Z£lT 1=1 



This completes the proof of Theorem 4.2.3 



□ 



The gist of the above theorem is the inequality 



(4.2.25) 



We will use it repeatedly in the proof of the following theorem. 

The following theorem also shows an order of \fa given that the state space of P is 



finite. It uses a similar approach as in the proof of Theorem 4.1.1 and is much lengthier 



than the previous result, but it gives a better coefficient of yfa in the sense that it does 
not depend on the size L. Therefore it gives some insight on the possibility to remove 
the finiteness condition. 

Theorem 4.2.4. Let {Fj}j & % + be an i.i.d. sequence under P from a finite set of distri- 
butions {Hi, H L }. Again assume for each t, 

H e (x) 1/2 dx + (l - H e (x)) 1/2 dx < oo. (4.2.26) 

-oo JO 

Then as a \ 0: 



^(1, a) < ji + 2a*y/a + o(y/a 



(4.2.27) 
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where //* = max^{/j(i/^)} and a* = y^max^{a 2 (H^)}. 

Proof. Again if we assume n(He) = fi* for all i = 1,...,L, \I/(l,a) will get no smaller. 
Without loss of generality, we will assume fi(He) = hereafter unless specified otherwise. 



In a similar way as we did in the proof of Theorem 4.1.1 we let r = r(a) be a positive 
integer-valued function such that r(a) /* oo and r^fa \ as a \ 0. Define the r x 1 
blocks as B r (x, y) = {{rx + i, y) : i — 0, 1, r — 1} for (x, y) G Z+. 

For every point z = G Z+, write X r (^) = ^ u6Br ( z ) X(v). The distribution 

function of X r (z) is denoted ), and the corresponding last-passage time function 

for X r (z) is then denoted as ^ r (x, y). Also, write for the convolution He* Hi* ...*Hi 
with He repeated r times. 

For j G Z + , let $ rj be the distribution function of the normal distribution 7V(0, rcrj), 
where <rj = Var(Fj). For each z = G Z+, let 3^(z) be a random variable with 



distribution $ rj -. Write 



r (x,y) = lim — E max V^y^z). 



n->co n TreU([nx\,[ny\) ■ 

2£7T 

Also, let $f be the distribution function of jV(0, rVf), where 14 is the variance of iifY 
We make an approximation first: 

Lemma 14. 

lim— a) - -tf*(l,ra)| = 0. (4.2.28) 

<4o Ja r 



Proof. For each £ = 0,1,..., L and each z = G Z+, we do the following cou- 

pling. Define {u(z) : z = G Z^_} be i.i.d. Uniform(0, 1) random variables. Set 

X r (z) = F r j~ l {u(z)), where F r j~ (u) = sup{x : F r j(x) < u). Similarly define Y r (z) = 
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$-j(u(z)). Also define X r/ {z) = F r j-\u(z))I {Fi=Hl} and Y r/ (z) = ^{u{z))I {Fj=Hl} . 

We first assume {Hi, H 2 , . . . , Hi} are uniformly bounded and directly quote the 
computation from Lemma 4.2 and Lemma 4.5 in [13J: 



lim-E max Vjrt-^))-^ 1 ^ 

Z£7T 

/oo 
|^(s)-$J(s)|3£fa 
-oo 

<2v/ra(l + ra)) J yjcr~s(l + \-^— \)-*ds 



<Cr±oi\/oi{\ + raj 
for proper constant C that are independent of r, £ and a. 



Then we apply (4.2.25): 



tf r (l,ro) - tf* r (l,ro) 
< Bm -E max (XJz) - Y r (z)) 



<lim-E max ^((^M*)) ~ (*f) _1 («(*))) 

n->oo n 7ren(n,LnroJ) z — ' v 

Z£7T 

L 

= Cr 3 (T^A/a(l + ra). 



£=0 



(4.2.29) 



ra— >oo n 7ren(n,LnraJ) 1 

Z£7T 

£: x 

< > lim -E max > (X r Az) — Y r£ (z)) , A 

-A^n^n weU{ \ nrai) A^\ r ^ i r ^ i) (4.2.30) 

1=0 Z&TT 



If we switch the two terms on the left hand side, we can repeat the calculation and 
obtain 

1 K 

— =|tf r (l,ro) - #*(l,ra)| < Vc^-*VT+ra, (4.2.31) 



k=0 

which goes to as a \ 0. 
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Lemma 4.4 in [TH] still applies here since X(z)'s are uniformly bounded, so we have 



lim— =|\|j(l,cO - -V r (l,ra)\ = 0. 
c40 J a r 



Therefore 



1 1 

lim-= *(l,a) - -** r (l,ra) 

<*4-0 Ja r 



0. 



(4.2.32) 



If we do not have uniform boundedness, Lemma 4.3 of [12] says that for any e > 0, 
we can find distribution functions Hi for each k with bounded support, and with mean 
and variance equal to those of Hi, and J^° oo \He(s) — He(s)\^ds < e. Hence ( |4.2.32 ) holds 
for a), the corresponding last-passage time function associated with {Hi, . . . , H L }. 



We repeat the argument leading to (4.2.29) and (4.2.30) and get 



|*(l,a) - *(l,a) 

L 



/oo ^ 



H,(s)\*ds 



<2LeWa(l + a). 



(4.2.33) 



We let e approach and it follows that (4.2.32) is also valid for {Hi}. The lemma is 
proved. □ 



With Lemma 4.2.32, we know that ^(1, a) has the same coefficient of the term y/a 



with -^$ r (l,ra). \&$ r (l, ra) is the last-passage constant of a model where all X(z) 
follow mean zero normal distributions. The following lemma looks at the role played by 
the variances of normal distributions. 



Lemma 15. Let X and Y be independent random variables and X ~ A/"(0, a 2 ), then 
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E(X V Y) is an increasing function of a. 
Proof. Firstly, we have 



/°° f y i r 2 r°° 1 r 2 



Note that 



we get 



i a: 2 f 00 1 , a: 2 



a: . exp( ^)c&r + / x . exp( -)dx = 0, 

v^a PV 2a 2 ; A. PV 2a 2 ; 



1 x 2 f°° 1 x 2 

f y 1 a; 2 H 1 a; 2 

/-„ ^ eXP <- 2^ )& - /_„ * ^ eXP< ~ ^ >dX 



r° i x 2 

= J (y-ax)-f=ex.p(-—)dx. 



The last line can be viewed as a function of a and its derivative is 



1 P , *\ , 

/ rr exp(— — )dx > 0, 

J-oo 2 



^ 2 

because when y < 0, it is easy to see J_ CT a;exp(— ^-)^ < 0; when y > 0, 
/j oo a;exp(-^)rfa; < /f^ x exp(-^)cte = 0. 

Therefore £(I V 7) is an increasing function of a. □ 
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With Lemma 15, we can run a similar argument with the one leading to (4.2.20). 



That is, we fix a zq G and claim that E max^- e n(n, [n«j ) 2~2 z <=n Yr{z) is maximized 
when Y r (zo) has the largest possible standard variance ra* 2 . Repeat this reasoning we 
see that an upper bound for ^/$ r (l, ra) is given if we let all sites z G Zi have the largest 



possible variance. So (1.0.3) can be applied here: as a \ 0, 



\P$ r (l, ra) < 2y/ra* y/ra + o(y/ra). 



From 4.2.32 as a \ 0, we have 



tf(l,ar) < 2a*v^ + o(v^). 



(4.2.34) 



Consider the case //* 7^ 0, the last result becomes 



a) < fi* + 2a*y/a + o(yfa) 



(4.2.35) 



□ 



Remark 16. This theorem did not remove the finiteness of the state space because the 



approximations (4.2.30) and (4.2.33) depend on the size L. A more accurate method of 



approximation is needed in order to lift this assumption. 
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